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I.  INTRODUCTION' 

One  of  the  damage  mechanisms  for  military  equipment  in  a  nuclear 
blast  environment  is  the  whole  body  movement  that  may  be  produced. 
Tumbling  or  overturning  of  vehicles  can  cause  significant  structural 
damage  or  render  them  inoperable  until  righted.  Even  if  overturning  does 
not  occur,  a  change  in  orientation  may  be  produced  that  can  put  radar  and 
communication  equipment  out  of  action  until  the  proper  orientation  is 
restored . 

Experimental  data  describing  the  movement  of  military  equipment 
during  exposure  to  blast  waves  simulating  those  from  nuclear  explosions  is 
needed  for  survivability  assessments  of  that  equipment.  In  addition,  such 
data  are  required  for  validating  computer  models  of  vehicle  whole-body 
movement  under  long-duration  blast  loading.  Validation  of  such  codes  is 
very  important  so  that  the  response  of  U.S.  and  foreign  vehicle  systems 
can  be  predicted  with  reliability.  It  is  not  practical  to  obtain  vehicle 
response  data  experimental!}  for  the  full  range  of  targets,  yields,  and 
surface  conditions  of  interest  to  the  Army. 

Experiments  have  been  conducted  previously  in  which  vehicles  were 
exposed  to  the  blast  waves  from  large  explosions.  The  primary 
instrumentation  used  on  such  tests  to  record  whole-body  movement  has  been 
high-speed  photography  by  cameras  mounted  on  towers  near  the  vehicles. 
Unfortunately  for  this  technique,  the  blast  wave  raises  dust  from  the 
ground  surface  which  quickly  obscures  the  view  of  the  target.  Only  at 
relatively  low  blast  overpressures  (less  than  70  kPa)  have  useful 
displacement  versus  time  data  been  obtained.  At  higher  overpressures,  the 
dust  obscures  the  targets  before  significant  movement  occurs.  Dust 
suppression  materials  have  been  applied  to  the  ground  surface,  but  so  far 
they  have  not  prevented  dust  obscuration.  Thus  no  useful  displacement¬ 
time  data  for  tracked  vehicles  has  been  obtained  using  photography,  and 
that  for  wheeled  vehicles  is  limited. 

Large  nuclear  blast  simulators  have  been  constructed  and  used  for 
blast  loading  of  vehicles.  However,  targets  in  such  simulators  may  be 
obscured  by  a  condensation  cloud  formed  by  water  vapor  and  possibly  carbon 
dioxide  in  the  driven  gas. 

A  means  other  than  photography  is  required  to  measure  the 
displacement  versus  time  of  a  vehicle  struck  by  a  blast  wave.  One  method 
of  interest  is  the  placement  of  a  self-contained  inertial  tracking  system 
on-board  the  vehicle.  Such  a  system  will  use  accelerometers  and  roll  rate 
sensors  whose  output  will  be  recorded  in  on-board  memories.  After  the 
test,  a  computer  program  will  be  used  to  derive  the  displacement  and 
rotation  versus  time  for  the  vehicle  from  the  recorded  data. 

The  objective  of  this  study  was  to  investigate  the  use  of  inertial 
reference  systems  for  measuring  blast-induced  displacements  of  military 
targets . 
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II.  PROCEDURE 


The  procedure  was  as  follows: 

U)  Review  the  literature  for  inertial  reference  systems  and  select 
an  appropriate  sensor  array. 

(1)  Design  and  construct  a  test  facility  to  translate  and  rotate  a 
sensor  array  at  approximately  the  rate  that  would  be  experienced 
during  the  overturning  of  a  military  vehicle. 

(3)  Test  and  evaluate  the  capability  of  a  sensor  array  of  off-the- 
shelf  components  to  describe  the  movement  produced  by  the  test 
facility. 

(4)  Develop  the  requirements  for  a  complete  inertial  reference 

system . 


III.  INERTIAL  REFERENCE  SYSTEMS 

The  primary  application  of  inertial  reference  systems  has  been  for 
navigation.  Gimballed  inertial  platform  systems  are  by  far  the  most 
common  type  of  inertial  system  in  military,  space,  and  commercial  use. 
Gyroscopes  and  accelerometers  on  the  stable  platform  are  inertially 
stabilized  from  the  body  motion,  and  the  stable  platform  physically 
represents  an  inertial  reference  frame.  By  double  integration  of  the 
specific-force  indications  from  the  accelerometers,  with  a  correction  for 
gravity,  position  determination  is  accomplished.  The  gimbals  isolate  the 
gyroscopes  from  the  possibly  large  rotation  rates  of  the  body  in  which  the 
system  is  placed.  The  state-of-the-art  of  gimballed  systems  is  such  that 
their  navigation  performance  can  approach  almost  error-free  instrument 
operation  to  a  point  where  uncertainties  in  the  knowledge  of  the  gravity 
field  become  the  dominant  source  of  navigational  error. 

Another  type  of  inertial  reference  system  is  the  strap-down  inertial 
system.  In  such  a  system,  the  sensors  are  mounted  directly  (or  perhaps 
with  vibration  isolators)  on  the  body.  The  inertial  sensor  outputs 
represent  specific  accelerations  and  angular  rates  coordinated  to  the  body 
axes.  To  function,  the  strap-down  systems  must  utilize  a  computer¬ 
generated  transformation  matrix  to  process  the  angular  rate  data  to  relate 
the  attitude  of  the  body  frame  axes  to  an  arbitrarily  specified  reference 
coordinate  system  with  a  known  relation  to  inertial  space.  This  attitude 
determination  is  then  used  in  the  transformation  of  the  acceleration 
information  to  the  reference  coordinate  system  for  derivation  of 
velocities  and  displacements.  A  detailed  discussion  of  strap-down  systems 
and  transformation  algorithms  is  given  in  Reference  1.1 


1  "Strap-Down  Inertial  Systems,"  AGARD  Lecture  Series  No.  95,  AGARP- 
l.S-93,  Advisory  Group  for  Aerospace  Research  and  Development,  North 
Atlantic  Treat'-  Organization,  7  rue  Ancelle,  92200  Neuilly  sur  Seine, 
trance,  Mav  1978.  ADA055778 


In  addition  to  the  added  computations,  the  inertial  sensors  require  a 
larger  dynamic  range  than  gimballed  systems.  Since  they  are  attached  to 
the  body,  they  are  subjected  to  the  entire  body  dynamic  environment,  and 
generally  will  not  perform  as  well  as  if  they  had  been  isolated  from  it  by 
gimbals.  However,  such  systems  have  been  pursued  for  navigation  purposes 
because  they  can  be  smaller,  easier  to  maintain,  and  of  less  cost. 

The  system  components  designed  for  navigational  purposes  are 
generally  unsuitable  for  use  in  the  blast  environment  because  of  their 
cost  and  susceptibility  to  shock.  The  applications  other  than  for 
navigation  found  in  the  literature  have  used  various  strap-down  systems. 

A  military  vehicle  struck  by  blast  may  move  with  six  degrees  of 
freedom  as  a  rigid  body.  The  translation  is  expressed  by  three 
independent  coordinates  which  specify  the  location  of  a  fixed  point  in  the 
rigid  body.  The  other  three  degrees  of  freedom  correspond  to  changes  in 
orientation  of  the  rigid  body  and  can  be  described  by  rotational 
parameters.  The  rotational  motion  can  be  described  in  terms  of  direction 
cosines,  quar ternions ,  or  Euler's  angles. 

Figure  1  shows  a  rigid  body  moving  in  space.  The  1,2,3  coordinate 
system  is  located  in  the  rigid  body  (body-centered  system)  and  moves  with 
it,  while  the  x,y,z  coordinate  system  is  the  fixed  system.  For  a  rigid 

body , 


r  =  R  +  uxj)  +  ajx  (w  x  j>)  (1) 

where  all  quantities  are  vectors  (underlining  denotes  a  vector  quantity), 

and 

r  =  acceleration  of  point  P  relative  to  origin  of  fixed  system 

K  =  acceleration  of  the  origin  of  the  body-centered  system 

=  angular  velocity  of  the  body-centered  system  relative  to  the  fixed 
system 

time  rate  of  change  of  angular  velocity 
[>  =  position  of  point  P  as  measured  in  the  body-centered  system. 

A  measured  acceleration  contains  a  component  due  to  gravity, 

a  =  r  +  £,  (2) 

where  a  =  acceleration  measured  by  an  accelerometer,  and 


g  =  acceleration  due  to  gravity. 


Then, 


a  =  R  + 


X£+(£X  (yx£)+£ 


(3) 


Transforming  vectors  from  the  body-centered  system  to  the  fixed 
system,  or  vice-versa,  requires  the  use  of  a  rotational  transformation 
matrix  and  translation  of  the  origin  of  the  coordinate  system. 2,3,4 


Once  the  acceleration  has  been  determined 
system,  velocity  is  calculated  as  follows: 


in  the  fixed  coordinate 


v  =  /(a  -  £)  dt  +  vc 


(4) 


Position  is  obtained  by  integrating  the  velocity: 


r  =  ,'v  dt  +  r. 


(5) 


The  initial  velocity  v0  and  initial  position  r0  must  be  known  in 
order  to  perform  the  integrations. 

A  variety  of  sensors  and  sensor  configurations  can  be  placed  in  the 
moving  rigid  body  to  form  a  strap-down  inertial  tracking  system.  The 
sensors  can  be  linear  accelerometers  only,  or  combinations  of  three  linear 
accelerometers  and  rate  sensors  and  angular  accelerometers.  The 
applications  found  in  the  literature  have  used  mainly  arrays  of  linear 
accelerometers . 

The  Naval  Biodvnamics  Laboratory  (NBDL)  has  pursued  a  program  to 
study  human  dynamic  response  to  various  stimuli  including  impact  and 


^  Mital,  N.K.  and  King  A. I.,  "Computation  of  Rigid-Body  Rotation  in 
Three-Dimensional  Space  From  Body-Fixed  Linear  Acceleration  Measurements," 
Transactions  of  the  ASME,  Journal  of  Applied  Mechanics,  Vol .  46,  December 
1979,  pp  925-930. 

3  McMaster,  D.K.,  "The  Development  of  an  Accelerometer  System  for 
Measuring  Pelvic  Motion  During  Walking,"  Ph.D.  Dissertation,  University  of 
Oxford,  1979.  AD-A088579 

^  Jeffers,  Jr.,  M.F.,  "Analytical  Methods  for  Determining  the  Motion 
of  a  Rigid  Body  Equipped  with  Internal  Motion-Sensing  Transducers," 
DTNSRDC-76-004 1 ,  David  W.  Taylor  Naval  Ship  R&D  Center,  Bethesda, 
Maryland,  October  1976. 
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'  In  tne  program  it  was  necessary  to  measure  the  acceleration 
:u:r.pone::t-.  In  a  typical  test  a  human  volunteer  was  instrumented 
! i rated  >>n  a  sled  that  impacted  a  barrier  that  produced  a 
:  deceit-ration  or  a  duration  of  0.25  to  0.30  seconds.  A 

inertial  tracking  system  for  installation  on  the  head  was 
•ft  used  a  planar  array  of  six  linear  accelerometers.  This 
re: v rence  s  'stem  was  iound  to  be  satisfactory  for  trajectories 
'...r  displacement  was  not  greater  than  a  quarter  of  a  revolution 
instrumented  subject  did  not  strike  obstructions  that  would 
ns.  rp  peaks  of  acceleration. 

requirement  for  minimal  mass  limited  the  choice  of  sensors,  and 
huted  a  careful  evaluation  of  two  types,  a  subminiature  piezo- 


•  inear  accelerometer1 


and  an  Endevco  angular  accelerometer. 


inear  accelerometers  were  usable,  but  only  if  they  were  carefully 
stringent  calibration  procedures  were  followed,  and  their 
■at  ire  maintained  constant  during  a  test.  It  was  found  that  the 
ir  a ...  ce  :  c  r  one  ter  also  should  be  maintained  at  the  same  temperature 
t-.-st  as  was  used  during  its  calibration,  and  should  be  protected 
■a:  id  temperature  changes. 


a.  nave  proceeded  to  an  inertial  reference  configuration  of 
accelerometers  in  a  plane.  Three  orthogonal  triads  are 

:he  vertices  of  a  triangle. 

linear  accelerometer  configuration  has  also  been  employed  by 
ling.-  The  head  rotations  of  anthropomorphic  dummies  and  of 
■e  measured  using  both  the  inertial  reference  system  and  high- 
;raphy.  In  the  tests,  the  dummy  or  cadaver  was  accelerated  at 
on  a  sled  and  impacted  a  device  that  produced  a  controlled 
The  frontal  impact  produced  a  peak  acceleration  of  about  18 
.si::. urn  rotation  of  the  head  was  about  120  degrees.  The 
■'  the  inertial  reference  system  results  to  those  from  high- 
:t  aphv  was  satisfactory.  Mital  and  King  concluded  that  they 
>ed  an  effective  and  reliable  method  to  compute  rigid-body 
three-dimensional  space  from  angular  accelerations  derived 
i.\cd,  nine  linear  accelerometer  reference  system.  The  study 
<■  advantage  of  computing  the  rotation  using  the  accelerometer 


■  ir.rj  '.villems,  <1 . ,  "An  Experimentally  Validated  3-D 
Package  tor  Application  in  Biodynamic  Research," 
I..:. eteent.ii  Stapp  Car  Crash  Conference,  San  Diego,  CA, 
,  pp  899-9  R). 


Detailed  Performance  Evaluation  of  Subminiature 
••i  I  ro  '■  t  e  i  s  , "  Proceedings  of  the  23  rd  ISA  international 
am ,  i->  May  1977,  pp  531-540. 

iet ailed  Evaluation  of  the  Endevco  Mode!  7302 
•  r  •  r  ,  "  !)i.-83p')()'j ,  Nava  1  Biodynamics  Laboratory,  New 


data  rather  than  measuring  them  from  film.  The  latter  procedure  was  slow 
and  very  laborious. 

McMaster^  developed  an  inertial  reference  system  using  a  planar  array 
of  eight  linear  accelerometers  (four  pair)  for  determining  the  three- 
dimensional  motion  of  the  pelvis  during  walking.  The  array  is  shown  in 
Figure  2.  He  derived  the  equations  for  the  eight  accelerometer,  four-pair 
configuration  and  the  method  of  solving  the  equations  for  angular 
velocity,  angular  position,  and  the  translational  acceleration,  velocity 
and  position.  He  found  it  necessary  to  construct  his  own  accelerometers. 
The  range  of  interest  was  +  5  g's,  and  the  maximum  digitization  rate 
required  was  200  Hz.  He  found  that  temperature  effects  on  the 
accelerometers  were  significant  and  had  to  be  minimized  by  special 
techniques.  However,  the  inertial  reference  system  provided  results  for 
walking  subjects  that  were  satisfactory  and  consistent  with  previous 
studies . 

Work  done  by  Jeffers^  at  the  David  W.  Taylor  Naval  Ship  R&D  Center 
seems  more  appropriate  for  the  problem  of  tracking  military  vehicles 
struck  by  blast.  The  Navy  interest  was  in  performing  model  launching 
studies  of  a  variety  of  torpedoes  and  missiles.  High-speed  photography 
was  an  available  technique  but  required  extreme  care  in  implementing,  and 
the  analysis  was  laborious  and  time  consuming.  Jeffers  derived  the 
appropriate  equations  for  two  sensor  configurations;  three  linear 
accelerometers  and  three  angular  rate  sensors;  and  five  linear 
accelerometers  and  one  angular  rate  sensor.  The  preferred  configuration 
was  three  accelerometers  and  three  rate  sensors.  The  formulation 
considered  body-oriented  axes  that  are  appropriate  for  a  military  vehicle. 
The  method  for  deriving  the  rotation  and  translation  parameters  from  the 
sensor  outputs  includes  accounting  for  the  initial  zeroing  of  signal  level 
in  the  gravity  field. 

Jeffers,  in  a  private  communication,  stated  that  limited  experiments 
were  performed  and  the  results  from  the  sensors  processed  using  his 
equations  agreed  well  with  results  from  photography.  The  angular  rate 
sensor  was  a  prototype  Humphrey  rate  sensor.  The  program  was 
discontinued,  and  Jeffers'  computer  code  for  processing  the  data  was  not 
documented  and  published. 

These  previous  studies  were  restricted  in  the  choice  of  sensors  by 
requirements  for  minimum  size  and  mass  to  permit  installation  of  the 
sensor  array  on  the  subjects  of  interest  and  to  reduce  the  influence  of 
the  sensor  array  on  the  response  of  the  subject.  Military  vehicles  are 
large  in  size  and  mass  compared  to  any  reasonable  sensor  array,  and  hence 
there  is  more  freedom  of  choice  for  sensors. 

The  analysis  of  data  from  a  six  accelerometer  array  was  found  to  be 
unstable. ^  The  eight  accelerometer,  four-pair  array  of  McMaster^  seems 
preferrable  to  the  nine  accelerometer  array  of  NBDL  or  of  Mital.^  it  uses 
one  less  accelerometer  and  does  not  place  more  than  one  accelerometer  at 
the  same  location.  However,  the  reported  experience  with  accelerometer 
arrays  indicates  a  need  for  very  accurate  gages  and  great  care  in 
calibration  and  in  temperature  control.  In  the  analysis  of  the  output 


ACCELEROMETER 

PAIRS: 


The  planar  eight  accelerometer,  four-pair  configuration  by 
McMaster  (Reference  3). 
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from  such  an  array,  the  difference  between  the  records  for  each  pair  of 
accelerometers  is  used  to  determine  the  rotational  components. 


Difficulties  may  occur  in  that  each  accelerometer  is  required  to 
sense  accurately  acceleration  due  to  translational  motion,  the  gravity 
field,  and  the  rotation  of  the  rigid  body.  The  magnitudes  involved  may  be 
incompatible  with  the  use  of  a  single  sensor  to  record  them. 

In  the  sensor  array  selected  by  Jeffers,^  consisting  of  three 
accelerometers  and  three  rate  sensors,  the  rotational  components  are 
determined  by  direct  measurement  of  angular  velocity  and  by 
differentiation  of  the  angular  velocity  to  obtain  angular  acceleration. 
For  maximum  reduction  of  error,  three  angular  accelerometers  can  be  added 
to  the  sensor  array  to  make  direct  measurements  of  angular  acceleration, 
thus  avoiding  the  necessity  for  differentiating  the  angular  velocity 
records.  However,  adding  these  accelerometers  increases  the  cost 
significantly  and  introduces  errors  from  the  accelerometers,  so  that  the 
reduction  in  error  achieved  may  not  be  worthwhile.  Although  the 
accelerometers  will  still  be  required  to  sense  accurately  the 
translational  and  gravitational  acceleration  components,  the  rotational 
components  will  be  derived  independently  from  the  rate  sensors  and 
possibly  angular  accelerometers.  Therefore  the  accelerometers  and  rate 
sensors  can  be  ranged  independently  of  each  other,  which  can  lead  to 
improved  accuracy. 

Based  on  this  review  of  previous  work,  the  selected  sensor  array  for 
an  inertial  tracking  system  for  military  vehicles  is  three  linear 
accelerometers  placed  in  a  triad,  if  possible,  and  aligned  with  the 
primary  axes  of  the  military  vehicle,  and  three  rate  sensors. 

The  potential  for  growth  of  error  seems  high  for  any  inertial  sensor 
array  under  blast  loading  conditions.  Before  beginning  the  construction 
of  an  inertial  tracking  system  capable  of  recording  the  full  six  degrees 
of  freedom  of  a  rigid  body,  it  seemed  necessary  to  show  that  tracking 
could  be  accomplished  with  satisfactory  accuracy  for  a  system  with  a 
minimum  of  sensors.  The  array  selected  for  evaluation  was  that  of  two 
linear  accelerometers  and  one  rate  sensor  placed  in  a  configuration  to 
track  motion  in  a  plane. 


IV.  DESIGN  AND  CONSTRUCTION  OF  TEST  FACILITY 


A .  Representative  Vehicle  Response 

One  military  vehicle  system  of  concern  for  blast  survivability  and 
hardening  studies  is  that  of  a  21/2  ton  truck  carrying  an  electrical 
equipment  shelter.  Such  a  system  has  been  exposed  to  blast  on  large-scale 
high  explosive  (HE)  tests  and  in  blast  simulators. 


A  calculation  of  the  overturning  of  such  a  system  by  blast  from  an  HE 
explosion  equivalent  to  a  one-kiloton  nuclear  explosion  was  available. 
The  code  used  was  B1.0M,  a  single  degree  -of-freedom  code  developed  at  the 
BRL  to  calculate  blast  overturning  of  vehicles.  The  blast  wave 


overpressure  was  68.9  kPa  (10  psi),  atmospheric  pressure  was  86.7  kPa 
(12.6  psi),  ana  atmospheric  temperature  was  12.0  degrees  Celsius.  The 
system  overturned  at  this  overpressure,  but  it  was  not  far  beyond  the 
limit  for  overturning. 

Figure  3  shows  the  angular  velocity  of  the  system.  It  has  an  initial 
peak  ol  80  degrees  per  second,  decreases  as  the  center  of  gravity  passes 
over  tiie  axis  of  rotation,  and  reaches  a  maximum  of  117  degrees  per  second 
at  impact  with  the  ground. 

Figure  4  shows  acceleration  verms  time  as  would  be  sensed  by  two 
accelerometers  located  at  the  center  of  gravity;  one  with  the  axis 
initially  vertical  in  the  plane  of  rotation,  and  the  other  initially 
horizontal  in  the  plane  of  rotation.  The  effect  of  the  gravity  field  on 
the  response  of  the  accelerometers  was  considered  in  generating  the 
curves . 

The  lower  set  of  curves  shows  an  initial  high  magnitude  pulse 
corresponding  to  the  diffraction  loading  period  followed  by  a  lower 
magnitude,  longer  duration  tail  that  corresponds  to  the  drag  loading 
phase.  After  0.19  seconds  the  system  is  rotating  under  the  influence  of 
gravity  alone,  and  the  time  until  impact  with  the  ground  is  long  compared 
to  the  positive  phases  of  acceleration.  The  magnitudes  of  the 
accelerations  during  this  period  are  small  compared  to  the  initial  peak 
magnitudes.  At  the  end  of  blast  loading,  the  system  has  rotated  13 
degrees.  Overturning  is  completed  at  1.24  seconds. 

The  upper  set  of  curves  in  Figure  4  show  that  the  initial  pulse  has  a 
duration  of  about  0.15  seconds.  The  calculation  assumed  a  perfectly  rigid 
body.  In  a  real  system  the  peak  accelerations  would  be  reduced  by  the 
flexure  of  the  structure  and  the  pulse  width  would  be  somewhat  broadened. 

To  approximate  overturning  of  the  truck-shelter  system  for  the  case 
described  above,  the  test  facility  should  provide  a  loading  that  produces 
an  initial  acceleration  pulse  of  about  10  g's  that  decreases  in  0.015- 
0.020  milliseconds  to  about  2  g's,  and  then  decays  to  zero  in  0.15  -  0.2 
seconds.  Specification  of  a  peak  acceleration  of  about  10  g's  instead  of 
the  higher  values  shown  in  Figure  4  is  an  attempt  to  account  for  the 
reduction  caused  by  the  flexure  of  a  real  system.  The  facility  should 
produce  overturning  in  about  1.2  seconds. 

H .  Concept  for  the  Test  Facility 

Figure  3  shows  the  concept  for  the  test  facility  to  produce  a 
rot  it  ion  of  90  degrees  about  an  axis  perpendicular  to  the  vertical  plane. 
The  rotating  component  is  a  quarter  wheel  constructed  to  minimize  its  mass 
and  moment  of  inertia  about  the  rotation  axis.  The  accelerometers  are 
placed  at  the  sensor  location  shown,  and  move  along  the  indicated 
tr.iject‘,1  as  the  quarter  wheel  rotates. 

: he  quarter  wheel  is  loaded  by  coil  springs  tied  to  it  by  a  cable, 
.'he  <  able  i  ;  wrapped  around  the  wheel  circumference  so  that  the  force 
moment  arr  with  respect  to  the  rotation  axis  remains  constant.  The 


Computed  angular  velocity  of  overturning  of  a  truck-shelter  combination  struck  by  a  68.9  kPa 
blast  wave  from  a  simulated  1  KT  nuclear  explosion. 


gure  .  Computed  accelerations  that  would  be  measured  by  two  accelerometers  placed  at  the  center 
an  overturning  truck-shelter  combination  struck  by  a  68.9  kPa  blast  wave  from  a  simulated 
1  KT  nuclear  explosion.  The  effect  of  gravity  on  the  accelerometers  is  included. 


Figure  5.  Concept  for  a  facility  to  simulate  overturning  of  a  military  vehicle  system  by  blast. 


quarter  wheel  is  held  in  place  by  the  restraint  cable.  To  execute  the 
test,  the  release  pin  is  pulled  free.  The  high  force  spring  unloads 
rapidly  and  simulates  blast  diffraction  loading,  and  the  low  force  spring 
simulates  drag  loading. 

C .  Modeling  of  Test  Facility 

To  provide  guidance  in  the  design  of  the  test  facility,  a  computer 
code  that  modeled  the  response  of  the  quarter  wheel  under  spring  loading 
was  written  and  utilized.  The  equation  of  motion  is  as  follows: 


where 


^axis 


t 

H 

K1 


spring  moment  -  weight  moment 

H[Kt  (Zy  -  Ho)  +  K2  (Z2  -  H0)]  -  W  Ry  cos  (6  +  0j) 

for  Zy  >  HO  and  Z2  >  HO,  and 

H[K2  (Z2  -  HO  )]  -  W  Rj  cos  (0  +  Qy)  for  Zy  <  H0 

and  Z2  >  H6 ,  and 

-W  cos  (0  +  yy)  for  Zy  <  H0  and  Z2  £  H0 , 


(6) 


d“e 


moment  of  inertia  about  axis 

angle  of  rotation  in  radians  from  the  initial  position  of 
the  quarter  wheel 

time 

constant  length  of  force  moment  arm 
spring  constant  for  high  force  spring 

maximum  extension  in  length  beyond  its  rest  position  of 
high  force  spring 

spring  constant  for  low  force  spring 

maximum  extension  in  length  beyond  its  rest  position  of 

low  force  spring 

radius  to  center  of  gravity  from  axis 
initial  angle  of  K . 


mass  of  wheel  and  sensors 


g 


acceleration  due  to  gravity 


W  =  weight,  mg 

The  mass  of  the  springs  is  small  compared  to  the  weight  of  the 
quarter  wheel,  and  is  neglected. 

Figure  6  shows  a  representation  of  the  model  of  the  quarter  wheel. 

To  compute  the  output  of  accelerometers  placed  on  the  quarter  wheel, 
the  components  of  acceleration  due  to  rotation  and  to  the  gravity  field 
must  be  considered.  Figure  7  shows  accelerometers  placed  at  the  radius  R2 
from  the  axis  of  rotation. 


At  the  radius  R2,  the  magnitudes  of  the  components  of  acceleration 
due  to  rotation  are 

ap  =  R2  0 

•2  (?) 

aN  =  8  2 

where 

ap  =  acceleration  magnitude  at  radius  R2  that  is  perpendicular  to  the 
radius  and  tangent  to  the  path  of  movement 

a^  =  acceleration  magnitude  at  radius  R2  that  is  along  radius  R2  and 
normal  to  the  path  of  movement 

R2  =  radius  of  accelerometers  from  axis  of  rotation 

0  =  angular  velocity 

••  d2e 

0  =  — r,  angular  acceleration 
dt 

If  the  accelerometers  are  oriented  initially  so  that  one  has  its 
sensitive  axis  vertical  and  the  other  horizontal,  they  will  sense  the 
components  of  acceleration  due  to  rotation  as  follows: 


aVR  =  aT  cos  ®2  “  aN  s*-n  ®2 
=  ap  sin  ©2  +  a^  cos  ©2 

where 


(8) 


ayR  =  acceleration  due  to  rotation  sensed  by  initially  vertical 
accelerometer 

=  acceleration  due  to  rotation  sensed  by  initially  horizontal 
accelerometer 

©2  =  initial  angle  of  radius  R2  with  respect  to  the  horizontal 
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Figure  7, 


Acceleration  components  at  the  location  of  the  accelerometers 
on  the  quarter  wheel.  Gravitational  acceleration  g  is  shown 
vertical  because  of  the  sign  convention  chosen. 


The  sign  convention  chosen  is  such  that  vertical  components  in  Figure 
7  that  would  cause  an  accelerometer  mass  to  move  down  are  called  positive, 
and  horizontal  components  that  would  cause  an  accelerometer  mass  to  move 
to  the  right  are  called  positive. 


If  the  vertical  accelerometer  is  set  to  read  zero  initially  in  the 
:ravity  field,  then  the  outputs  of  the  accelerometers  due  to  rotation  in 
he  gravity  field  are  computed  as  follows: 


and 


aV  =  aVR 

-  8 

(i 

-  cos  9) 

aH  =  aHR 

-  8 

sin 

1  6, 

av  =  aRT 

cos 

e2 

-  aRN  sin  &2  ~ 

aH  =  aRT 

sin 

92 

+  aR;>j  cos  02  - 

ay  =  83^ 

cos 

62 

-  R2  9^  sin  02 

aH  =  R?* 

sin 

e2 

+  R2  9^  cos  02 

(9) 


whe  re 


ay  =  acceleration  sensed  by  the  accelerometer  that  is  initially 
vertical 


aR  =  acceleration  sensed  by  the  accelerometer  that  is  initially 
horizontal 


Figure  5  shows  the  configuration  selected  for  the  quarter  wheel,  with 
aluminum  sheet  and  aluminum  channel  used  for  construction.  The  maximum 
radius  of  the  wheel  as  finally  designed  was  1.22  metres,  with  a  computed 
mass  (with  sensors)  of  18.1  kilograms  and  a  moment  of  inertia  about  the 
axis  of  15.25  kilogram  (metres)^. 


Two  cases  were  finally  selected  to  attempt  to  achieve  in  the  test 
facility,  Case  1  with  peak  accelerations  of  10  g’s,  and  Case  2  with  peak 

accelerations  of  5  g's. 

For  Case  1  and  Case  2,  H  =  1.194  metres,  =  0.862  metres,  R2  = 

!  .'*37  metres,  and  =  ©2  =  45  degrees. 

For  Case  1 ,  Kj  =  73550  Newtons/metre,  Zj  =  0.0127  metres,  K2  =  1576 
h'ewtons/metre,  and  Z2  =  0.254  metres. 

For  Case  2,  Kj  =  38530  Newtons/metre,  7,[  =  0.0127  metres,  K2  =  437.8 

Newtons/metre,  and  Z 2  ~  0.508  metres. 


Figure  8  shows  the  predicted  acceleration  versus  time  waveforms  for 
Cases  1  and  2  for  initially  vertical  and  horizontal  accelerometers  located 
at  radius  ks  at  an  initial  angle  of  45  degrees.  The  primary  waveform 


ACCELEROMETER  INITIALLY 
VERTICAL 


Predicted  accelerometer  waveforms  for  Cases  1  and  2  for  the  test  facility 


characteristics  of  a  hit'll  magnitude  initial  pulse  of  the  desired  magnitude 
and  duration  : o l lowed  by  a  low  magnitude  long  duration  tail  are  predicted 
for  Case  !.  Case  1  predicts  waveforms  i'or  a  higher  yield  simulation, 
where  overturning  occurs  at  lower  initiil  g  levels  and  a  longer  duration 

wave . 


!).  Construction  and  Assemble 

i he  test  facility  was  designed  for  construction  from  sheet,  channel 
and  angle  steel  and  aluminum.  The  radius  of  the  quarter  wheel  is 
nominally  0.12  metres  (4  feet).  Figure  9  shows  the  facility  as 
constructed  and  assembled.  The  basic  frame  is  made  of  steel  channel.  The 
quarter  wheel  is  held  in  the  vertical  plane  by  guy  cables  on  each  side 
that  are  attached  to  the  rotating  axle.  The  springs  that  apply  the  load 
are  attached  to  long  threaded  rods  that  pass  through  a  threaded  plate  at 
the  end  of  the  frame.  The  springs  are  extended  to  apply  the  loading  by 
unscrewing  the  threaded  rods.  The  wheel  is  held  in  place  by  a  restraint 
cable  that  in  turn  is  held  by  a  quick-release  pin. 

Figure  10  shows  the  quarter  wheel  in  its  initial  position.  The 
restraint  cable  is  shown,  as  well  as  the  guy  cables  and  their  attachment 
to  the  rotating  axle. 

Figure  11  shows  the  quarter  wheel  in  the  overturned  position.  After 
rotating  90  degrees  the  wheel  strikes  packing  foam.  The  magazines  seen  in 
the  figure  are  a  convenient  means  to  adjust  the  heights  of  the  initial 
support  and  final  impact  area. 

Figure  12  shows  a  rear  view  of  the  facility.  The  loading  cable  and 
its  attachment  on  the  wheel  can  be  seen.  The  wrap-around  length  of  cable 
is  0.79  metres. 

Figure  1  l  shows  the  pin  that  anchors  the  restraint  cable.  When  the 
springs  have  Seen  stretched  so  that  the  desired  force  is  acting  upon  the 
quarter  wheel,  this  pin  is  pulled  tree  and  wheel  rotation  begins. 


V.  INSTRUMENTATION' 


A.  Selection  of  Sensors 


For  evaluation  of  the  capability  of  an  inertial  reference  system  for 
tracking,  movement  in  a  vertical  plane,  two  linear  accelerometers  and  one 
gage  to  measui e  angular  velocity  were  required  for  mounting  on  the  quarter 
wheel.  The  sensors  found  for  measuring  angular  velocity  were  rate 
gyroscopes,  integrating  angular  accelerometers,  and  a  rate  sensor 
manufactured  b\  Humphrey,  Inc. 


Mil-  iiumphrev  rat 
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r.i  t  e  .  I  >e  f  1  •  •  c  t  i  on  ■  >  t 
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flow  over  one  of  two  fine 
of  the  flow  channel.  The 
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Figure  9.  The  overturning  test  facility. 


Figure  10.  The  quarter  wheel  in  the  initial  position. 


f i ow  past  a  wire  cools  the  wire  and  lowers  its  resistance,  unbalancing  a 
bridge  circuit  containing  both  wires.  The  change  in  the  bridge  is  related 
to  the  angular  rate  of  rotation.  This  sensor  offered  the  best  combination 
of  features  such  as  higher  natural  frequency,  damping,  accuracy,  and 
resistance  to  shock  and  acceleration.  The  specifications  for  the  Humphrey 
sensor  are  shown  in  Table  1. 

The  accelerometers  selected  were  closed-loop  force  balance 
transducers.  These  accelerometers  typically  have  errors  an  order  of 
magnitude  or  more  less  than  conventional  accelerometers.  The  particular 
model  purchased  was  the  Schaevitz  LSB  linear  accelerometer.  In  operation, 
a  pendulous  mass  in  the  accelerometer  develops  a  torque  proportional  to 
the  product  of  its  mass  unbalance  and  the  applied  acceleration.  The 
motion  of  the  mass  is  detected  by  a  position  sensor.  The  output  of  the 
position  sensor  is  applied  to  a  servo  amplifier  driving  a  torque  motor, 
and  a  torque  equal  and  opposite  to  the  acceleration  torque  is  applied  so 
that  only  minute  motion  of  the  mass  occurs.  The  current  through  the 
torque  motor  is  proportional  to  the  input  acceleration,  and  is  used  to 
produce  the  output  voltage.  Table  2  lists  the  specifications  for  the 
Schaevitz  accelerometers  used  in  this  study. 

Standard  strain  gage  accelerometers  manufactured  by  the  Endevco 
Corporation  were  also  available.  Table  3  lists  the  specifications  for  the 
model  used,  which  employs  piezo-resistive  elements  in  a  bridge  circuit. 

Two  other  types  of  sensors  were  installed  on  the  test  facility.  A 
precision  potentiometer  was  placed  on  the  rotating  axle  of  the  wheel,  and 
a  load  cell  was  placed  in-between  the  springs  loading  the  wheel  and  the 
cable  connecting  them  to  the  wheel. 

The  potentiometer  was  placed  in  an  electrical  circuit  so  that  a 
voltage  output  was  generated  proportional  to  the  rotation  of  the  axle,  and 
the  output  was  calibrated  in  terms  of  rotation  angle.  The  signal  from 
this  sensor  was  intended  to  provide  position  in  time  and  space  of  the 
accelerometers  for  reference  and  comparison  with  that  derived  from  the 
accelerometers  and  rate  sensor  records. 

The  load  cell  was  used  to  indicate  the  force  acting  to  turn  the 
quarter  wheel.  It  enabled  tests  to  be  set  up  for  the  acceleration  levels 
of  interest  and  to  repeat  tests  with  the  same  loading.  Its  output 
provided  an  unloading  versus  time  waveform.  Table  4  lists  the 
specifications  for  the  Kulite  load  cell  used.  The  cell  is  a  flat  load 
tell  employing  a  stainless  steel  diaphragm  construction.  The  transduction 
element  is  a  four-arm  Wheatstone  semiconductor  bridge,  internally 
i ompensuted  for  thermal  changes. 

e alibration  of  Sensors 

The  accelerometers  were  statically  calibrated  with  the  aid  of  a  Gyrex 
pin  tabie  located  at  the  BRL.  The  spin  table  was  capable  of  applying  to 
test  specimens  a  centrifugal  force  of  known  magnitude.  The  acceleration 
applied  to  the  accelerometers  was  calculated  from  the  following  relation: 


TABLE  1.  SPECIFICATIONS  OF  THE  HUMPHREY,  INC.,  RATE  SENSOR 


Model  Number: 

Serial  Number: 

Range : 

Natural  Frequency: 
Damping: 

Input  Voltage: 

Input  Current: 
Calibration  at  12  volts. 


RT03-0136-1 

101 

+  200  degrees/second 
25  Hz  minimum 
0.7  +  0.2 

+12+2  volts  D.C. 

200  milliamperes  maximum 
165  milliamperes: 


RATE  SENSOR  OUTPUT  DEVIATION 


+200 

deg/sec 

+2.524 

volts 

+.000 

volts 

+  160 

deg/sec 

+2.015 

volts 

+  .006 

volts 

+  120 

deg/sec 

+1.506 

volts 

+.012 

volts 

+80 

deg/ sec 

+1.004 

volts 

+.012 

volts 

+40 

deg/ sec 

+0.501 

volts 

+  .012 

volts 

+0 

deg/sec 

+0.001 

volts 

+  .009 

volts 

-40 

deg/sec 

-0.498 

volts 

+.006 

volts 

-80 

deg/sec 

-1.001 

volts 

+.006 

volts 

-120 

deg/ sec 

-1.498 

volts 

+.000 

volts 

-160 

deg / sec 

-1.998 

volts 

-.003 

volts 

-200 

deg/sec 

-2.503 

volts 

-.000 

volts 

Linearity  (Calculated) ; 
Scale  Factor  (Calculated): 
Noise : 

Null  Signal: 

Acceleration : 

Shock: 

Temperature : 

Mass : 


0.237  percent  of  full  scale 
0.0126  volts/degree/second 
10  millivolts  RMS 
+  50  millivolts 
100  g's,  any  axis 
100  g's,  11+1  millisecond 
-45  degrees  C  to  +77  degrees  C 
0.31  kilograms  nominal 


45 


TABLE  2.  SPECIFICATIONS  FOR  THE  SCHAEVITZ  LINEAR 
ACCELEROMETER  AT  20  DEGREES  CELSIUS 


Model : 

Range : 

Nominal  Natural  Frequency: 

Damping  Ratio: 

Nominal  Output  Impedance: 

Input  Voltage: 

Input  Current : 

Full-Range  Open-Circuit 
Output  Voltage: 

Linearity : 

Hysteresis : 

Resolution: 

Static  Cross  Axis 
Sensitivity: 

Bias : 

Sensitive  Axis  to  Case 
Alignment : 

Noise  Output: 

Operating  Temperature: 

Thermal  Coefficient  of 
Sensitivity: 

Shock  Survival: 

Mass : 


LSBC-10  Linear  Accelerometer 

+  10  g's 

140  Hz 

0.6  nominal 

2500  ohms 

+  15  volts  D.C.  +  10  percent  nominal 
25  milliamperes  D.C.  maximum 

+  5.0  volts  D.C.  +  1  percent 
+0.1  percent  of  full  scale* 

0.02  percent  of  full  scale  maximum 
0.001  percent  of  full  scale  maximum 

+  0.002  g  per  g 

0.1  percent  of  full  scale  maximum 
+  1  degree 

5  millivolts  RMS  maximum 
-55  degrees  C  to  +95  degrees  C 

0.02  percent  per  degree  C  maximum 
100  g's  for  11  milliseconds 
0.085  kilograms 


*Full  scale  is  defined  as  "minus  full  range  to  plus  full  range" 


TABLE  3.  SPECIFICATIONS  FOR  THE  ENDEVCO  ACCELEROMETERS 
AT  24  DEGREES  CELSIUS 


Model : 


Range : 


Sensitivity  at  10  volts 
excitation,  reference  100  Hz: 


Non-Linearity  &  Hysteresis: 
(Percent  of  reading,  maximum, 
to  full  range) 


Frequency  Response: 

(  +  5  percent  maximum,  reference 
100  Hz) 


Mounted  Resonant  Frequency: 


Damping  Ratio: 


Transverse  Sensitivity  (Max): 


Thermal  Sensitivity  Shift: 
(Reference  24  degrees  C) 


Zero  Measured  Output  (Max) : 


Thermal  Zero  Shift  (Max): 
(-18  to  +93  degrees  Celsius) 


Excitation: 


Input  Resistance: 


Output  Resistance: 


Acceleration  Limits 
(in  any  direction): 


Temperature  (Operating) : 


2262-25 


±  25  g's 


20  millivolts  per  g  typical, 
16  millivolts  per  g  minimum 


+  1  percent  (above  25  g's,  +  3  percent 
of  reading,  maximum,  to  50  g's) 


0  to  650  Hz 


2500  Hz 


+  0.15 


-  0.10 


3  percent 


Percent 


Degrees  C 


+  25  millivolts 


+  20  millivolts 


10.0  volts  D.C. 


1800  ohms 


1400  ohms 


0.028  kilograms 


250  g's  static  and  sinusoidal, 
2000  g's  shock 


■18  to  93  degrees  Celsius 


TABLE  4.  SPECIFICATIONS  FOR  THE  LOAD  CELL  OF 
KULITE  SEMICONDUCTOR  PRODUCTS,  INC. 


S# 

88 


Model : 

Range : 

Diaphragm  Thickness: 

Nominal  Deflection: 

Natural  Frequency: 

Overload : 

Operational  Mode: 

Rated  Electrical  Excitation: 
Max  Electrical  Excitation: 
Input  Resistance: 

Full  Scale  Output: 

Output  Resistance: 

Residual  Unbalance: 

N'on-L  inear  ity : 

Hysteresis : 

Repeatability: 

Resolution : 

Operating  Temperature  Range: 

Compensated  Temperature 
Range : 

Thermal  Effect  on  Zero: 

Thermal  Effect  on 
Sensitivity: 

Diameter : 

!'h it  kness  : 

m  t  i  n : 


TC-2000 

+  2224  Newtons  (500  pounds) 

1.65  millimeters  (0.065  inches) 
0,0076  millimeters  (0.0003  inches) 

12  kHz 

150  percent  of  rated  capacity 
Tension  and  compression 
10  volts  DC/AC  (RMS) 

15  volts  DC/AC  (RMS) 

3000  ohms 

75  millivolts  nominal 
1000  ohms 
+  2  percent 

0.5  percent  of  full  scale  output 
0.1  percent  of  full  scale  output 
0.1  percent  of  full  scale  output 
Infinite 

-40  degrees  C  to  120  degrees  C 

25  degrees  C  to  80  degrees  C 

+  3  millivolts  per  55  degrees  Celsius 

+  1  percent  of  full  scale  output/ 

55  degrees  C 

51  millimeters 

27  mill imeters 

Center  Stud 
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a  =  2.840  x  10“5  N2  r 


(10) 


where 


a  =  acceleration  in  g's 

N  =  speed  in  revolutions  per  minute 

r  =  radius  to  the  seismic  mass  in  inches 

Each  accelerometer  in  turn  was  subjected  to  accelerations  in  5  g 
increments  up  to  its  maximum  rated  range.  The  calibration  data  obtained 
agreed  with  those  supplied  by  the  manufacturer. 

The  calibration  data  supplied  by  the  manufacturer  was  used  for  the 
rate  sensor. 

The  potentiometer  for  measuring  the  rotation  angle  was  calibrated  by 
recording  the  change  in  output  voltage  for  a  rotation  of  the  quarter  wheel 
through  90  degrees. 

C .  Installation  of  Gages  on  the  Test  Facility 

Figure  14  shows  an  Endevco  accelerometer  and  the  Humphrey  rate  sensor 
mounted  on  the  quarter  wheel.  The  accelerometer  was  mounted  on  a  Masonite 
fiberboard  panel  6.35  millimeters  thick.  The  second  accelerometer  was 
mounted  on  the  opposite  side  of  the  board  so  that  the  sensitive  masses  jf 
the  two  accelerometers  were  at  the  same  radius  and  angle.  Several 
different  mounting  configurations  were  used  in  attempting  to  minimize 
mechanical  vibrations  at  the  accelerometers.  The  final  test  runs  were 
made  with  the  accelerometers  mounted  in  an  aluminum  bracket,  with  a  layer 
of  lead  and  a  layer  of  neoprene  rubber  between  the  aluminum  bracket  and 
the  fiberboard.  The  neoprene  rubber  was  placed  against  the  fiberboard. 

The  rate  sensor  was  mounted  directly  to  the  aluminum  channel. 

Figure  15  shows  the  potentiometer  that  was  coupled  to  the  end  of  the 
axle  to  measure  rotation  angle. 

Figure  16  shows  the  installation  of  the  load  cell  between  the  loading 
springs  and  the  cable  that  was  attached  to  the  quarter  wheel. 

D .  Data  Acquisition  and  Reduction  System 

The  block  diagram  of  a  typical  channel  of  the  data  acquisition  and 
reduction  system  is  shown  in  Figure  17.  The  B&F  signal  conditioners 
provided  excitation  voltage  and  balancing  for  the  Endevco  accelerometers. 
They  also  provided  operating  power  for  the  Schaevitz  accelerometers.  The 
Kulite  load  cell,  the  Humphrey  rate  sensor,  and  the  potentiometer  on  the 
axle  were  powered  by  separate  power  supply  modules.  The  signals  Iron  all 
the  transducers  were  amplified  with  a  Tektronix  AM  502  amplifier  be! .  re 
being  recorded. 
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The  Humphrey  rate  sensor  and  an  Endevco  accelerometer  mounted  on  the  quarte 


sion  poten 


Figure  16.  The  load  cell  placed  to  read  the  force  applied  by  the  springs 
to  the  cable  attached  to  the  quarter  wheel. 


Recorder 


A  Honeywell  Model  10 L  FM  cape  recorder  was  used  to  record  the 
signals.  It  was  configured  to  operate  in  the  intermediate  band  mode  at  a 
tape  recording  speed  of  0.381  metres  (15  inches)  per  second.  This 
configuration  provided  the  optimum  signal-to-noise  ratio  for  recording 
data . 

Calibration  steps  were  recorded  on  the  tape  prior  to  each 
experimental  session.  The  voltage  substitution  method  was  used.  In  this 
method,  a  precise  voltage  representing  the  known  value  of  the  parameter  is 
recorded  on  the  tape  for  each  of  the  data  channels. 

The  data  recorded  on  tape  were  played  back  through  a  Tektronix  AM  502 
amplifier  into  a  Pacific  Instruments  Transient  Digitizer  at  the  same  tape 
speed  as  was  used  for  recording.  The  data  were  examined  initially 
unfiltered.  For  input  into  a  computer  code  to  calculate  the  trajectory  of 
the  sensor  array,  the  data  were  passed  through  a  low-pass  filter  in  the 
amplifier  set  at  100  Hz,  and  redigitized.  These  data  were  transferred  to 
the  IBM  AT  computer  for  plotting,  processing,  and  permanent  storage  on  a 
magnetic  disk. 


VI.  PROGRAM  TO  COMPUTE  TRAJECTORY  OF  SENSORS 

A  program  was  written  to  compute  the  vertical  and  horizontal 
displacement  versus  time  of  the  sensor  array.  The  program  listing  is 
contained  in  Appendix  A. 

The  program  assumes  that  the  accelerometers  are  initially  horizontal 
and  vertical,  that  they  are  coincident,  and  that  the  vertical 
accelerometer  initially  is  biased  so  that  its  output  is  zero  in  the 
gravitv  field.  Figure  18  shows  the  geometry  involved. 

In  the  program,  the  angular  velocity  from  the  rate  sensor  is 
integrated  to  obtain  the  angle  of  rotation.  This  angle  is  used  to  correct 
the  output  from  the  accelerometers  for  the  effect  of  gravity.  It  is  used 
also  to  calculate  the  horizontal  X  and  vertical  Y  components  of 
acceleration  from  the  corrected  accelerometer  records.  The  X  and  Y 
components  of  acceleration  are  integrated  to  obtain  X  and  Y  velocities, 
and  integrated  again  to  obtain  X  and  Y  displacements.  The  calculation  is 
repeated,  using  the  rotation  angle  as  determined  from  the  potentiometer  on 
the  axle,  to  obtain  displacements  without  integrating  the  rate  sensor 
record  to  obtain  the  angular  velocity. 

Integration  is  accomplished  at  each  time  step  utilizing  the 
trapezoidal  rule.  The  only  control  of  the  accuracy  of  the  integration  is 
by  variation  of  the  time  step  value. 

The  program  begins  bv  reading  information  describing  the  test,  the 
gages  used,  and  spring  data.  Then  it  reads  data  required  for  the 
computation.  These  are  the  radius  from  the  rotation  axis  to  the  sensitive 
masses  of  the  accelerometers,  and  the  horizontal  distance  from  the  axis  to 
the  vertical  line  through  the  sensitive  masses  of  the  accelerometers. 


From  these  values  the  initial  angle  of  the  radius  to  the  accelerometers  is 
computed . 

The  sensor  records  fed  to  the  program  must  be  digitized  at  the  same 
rate  so  that  the  time  between  points  is  constant  for  all  records.  The 
ordinates  of  the  records  are  read  by  the  program,  and  the  time  for  each 
point  is  calculated  from  the  input  value  of  the  time  step.  Another  input 
is  the  number  of  the  time  step  at  which  the  computation  is  to  stop. 

The  gravitational  constant  is  defined  in  the  program  as  9.800  m/s^  at 
Aberdeen,  Maryland. 

To  begin  the  computation,  the  initial  values  of  the  parameters  are 
defined.  All  velocities  and  displacements  are  set  to  zero.  Starting 
values  for  the  accelerometers  are  computed  as  follows: 


aHgl  "  aHl 

(ID 

aVgl  =  aVl 

AX1  =  (aVgl  sin0 l  +  a^gi  cos©})  g  =  aHgig 
AY1  =  (aVgl  cos0L  -  aHgi  sin0 i)  g  =  aVgig 

where 

01  =  initial  value  of  rotation  angle  =  0 

g  =  acceleration  due  to  gravity 
am  =  starting  acceleration  for  the  horizontal  accelerometer 

ayi  =  starting  acceleration  for  the  vertical  accelerometer 

aHgl  =  aHl  corrected  for  gravity 

aVgl  =  aVl  corrected  for  gravity 

Ay i  =  starting  horizontal  (X)  component  of  acceleration  at 
accelerometers 

Ay i  =  starting  vertical  (Y)  component  of  acceleration  at 
acce lerometers 

Accelerations  described  by  lower  case  letters  are  in  terms  of  g's. 


Next  the  angular  velocity  is  integrated: 


start  time 


C1  = 

t2  =  second  time  value 

0rs1  =  initial  angular  velocity  =  0 
0RS2  =  angular  velocity  at  time  t2 

0RS1  =  initial  angle  of  rotation  =  0 

3rS2  =  integrated  value  for  angle  of  rotation  at  time  t2 

The  accelerometer  records  are  corrected  for  the  gravity  field  using 
the  angle  9rs2" 

aHg2  =  aH2  +  sin0RS2 

(13) 

aVg2  =  aV2  +  1  ”  cos6RS2 

where 

aR2  =  acceleration  at  time  t2  for  the  horizontal  accelerometer 

ay2  =  acceleration  at  time  t2  for  the  vertical  accelerometer 

aHg2  =  aH2  corrected  for  gravity 

aVg2  “  aV2  corrected  for  gravity 

The  X  and  Y  components  of  acceleration  at  time  t2  are  computed  from 
the  accelerometer  records  using  the  angle  ®rs2: 

aX2  =  (avg2  sin0RS2  +  aHg2  cos0RS2>  8 

(14) 

aY2  =  <aVg2  cos0RS2  _  aHg2  sin0RS2>  8 

where 

Ay 2  =  horizontal  (X)  component  of  acceleration  at  accelerometer  at 

time  t2 

Ay2  *  vertical  (Y)  component  of  acceleration  at  accelerometer  at  time 

t2. 

These  accelerations  are  integrated  to  obtain  X  and  Y  velocities: 

UX2  =  (1/2)  (aX2  +  AXl)  (c2  ~  cl)  +  UX1 

(  1  ;> ' 

UY2  =  (1/2)  (Ay 2  +  Ayj)  (c2  “  cl)  +  UY 1 

wliere 
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X;  *  starting  horizontal  displacement  of  the  accelerometers. 

The  value  is  zero. 

\>  -  horizontal  displacement  of  the  accelerometers  at  time  t£ 

V-  -•  starting  ve>-t  ical  displacement  of  the  accelerometers. 

The  .  t 1 ue  is  zero. 

-••*rt  :  c  <1  displacement  of  the  accelerometers  at  time  to 

i  i  integration  process  is  repeated  bv  the  program  using  the  angle 
•  >)  r  'tatl.T.  from  the  potent  iometer  record  instead  of  the  angle  obtained  by 
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These  results  for  the  first  integration  cycle  are  printed.  The 
computed  values  are  set  by  the  program  as  the  initial  values  of  the  next 
cycle,  and  the  sequence  cf  computations  is  repeated  until  either  the 
number  of  the  time  step  at  which  the  computation  is  to  stop  has  been 
reached  or  the  last  data  points  have  been  reached,  whichever  occurs  first. 

The  program  provides  three  descriptions  of  the  displacement  versus 
time  of  the  accelerometers.  The  most  accurate  is  that  derived  using  the 
potentiometer  record  and  the  radius  from  the  axis  to  the  accelerometers. 
The  second  and  that  of  primary  interest  is  derived  from  the  accelerometer 
and  rate  sensor  records.  The  third,  derived  using  the  accelerometer  and 
potentiometer  records,  provides  results  where  the  dominant  errors  should 
be  from  the  accelerometer  records. 


VII.  TEST  PROGRAM 

The  test  facility  provided  capability  of  movement  of  about  two  metres 
with  a  rotation  of  90  degrees  for  an  inertial  sensor  array.  The  sensor 
array  for  the  tests  consisted  of  two  accelerometers  and  the  Humphrey  rate 
sensor.  Two  types  of  accelerometers  were  used.  One  was  the  Endevco 
accelerometer,  a  piezo-resistive  strain  gage  type  with  a  natural  frequency 
of  2500  Hz.  The  other  was  the  Schaevitz  servo  accelerometer,  a  closed- 
loop  type  with  a  natural  frequency  of  about  140  Hz.  These  gages  are 
described  in  more  detail  in  Section  V.A. 

Two  loading  functions  were  used  in  the  program,  one  producing  a  peak 
acceleration  of  about  10  g’s,  and  the  other  about  5  g's.  Two  different 
sets  of  springs  were  used  to  generate  the  loading  on  the  quarter  wheel. 
Both  loading  functions  were  used  in  testing  sensor  arrays  containing  each 
type  of  accelerometer. 

A.  Test  Procedure 

Electrical  power  to  the  sensors  was  turned  on  at  least  five  minutes 
prior  to  a  test  for  sensor  warming  and  stabilization.  The  maximum  warm-up 
time  listed  by  the  manufacturer  for  any  gage  was  two  minutes. 

Loading  on  the  quarter  wheel  was  provided  by  a  long  spring  and  a 
short  spring  connected  in  parallel  in  the  cable  line.  First  the  long 
spring  was  connected  and  the  cable  centered  on  the  circumference  of  the 
wheel.  Then  the  restraint  cable  was  connected.  The  long  spring  was 
stretched  to  its  maximum  extension  by  unscrewing  the  long  threaded  rod  to 
which  it  was  attached.  The  force  in  the  cable  was  read  from  a  meter 
monitoring  the  output  of  the  load  cell.  The  spring  was  stretched  until 
the  force  in  the  cable  became  that  appropriate  for  the  loading  function 
that  the  spring  set  was  intended  to  produce.  This  was  either  400  Newtons 
(90  pounds)  or  222  Newtons  (50  pounds). 

The  short,  strong  spring  then  was  attached  to  the  cable  pulling  on 
the  wheel  and  stretched  by  unscrewing  the  threaded  rod  to  which  it  was 
attach ed.  This  spring  was  stretched  until  the  load  cell  monitor  showed 
the  desired  peak  force  in  the  cables.  The  value  was  either  1334  Newtons 
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t,300  pounds),  if  the  long  spring  loading  was  400  Newtons  (90  pounds);  or 
712  Newtons  (160  pounds),  if  the  long  spring  loading  was  222  Newtons  (50 
pounds ) . 

With  the  wheel  under  full  loading,  a  sensitive  bubble  level  was  used 
with  turnbuckle  adjustments  in  the  restraint  cable  to  set  the 
accelerometers  horizontal  and  vertical.  A  plumbob  was  used  to  define  a 
vertical  line  through  the  sensitive  masses  of  the  accelerometers.  The 
horizontal  distance  from  the  rotation  axis  to  this  vertical  line  was  read. 
The  radius  from  the  axis  to  the  accelerometers  was  measured.  These  two 
quantities  were  used  to  compute  the  initial  angle  of  the  radius  to  the 
accelerometers . 

Calibration  data  for  the  sensors  were  fed  into  the  tape  recorder. 
The  potentiometer  was  calibrated  by  rotating  the  wheel  90  degrees. 

For  recording  a  test,  the  tape  recorder  was  activated.  When  the  tape 
was  running  at  the  proper  speed,  a  marker  signal  was  placed  on  the 
potentiometer  channel  by  a  manually  operated  switch.  The  pin  holding  the 
restraint  cable  was  yanked  free,  the  springs  contracted  and  rotated  the 
quarter  wheel. 

The  trigger  signal  on  the  potentiometer  channel  was  used  as  zero  time 
for  all  channels.  The  records  were  digitized  as  recorded,  and  plotted  for 
evaluation.  Then  the  records  were  played  back  through  a  100  Hz  low  pass 
filter  and  digitized  at  about  0.5  millisecond  intervals.  The  baseline  for 
a  waveform  was  set  using  an  average  of  points  taken  shortly  before  the 
initial  rise  of  the  waveform.  These  data  were  fed  into  the  program 
described  in  Section  VI  and  the  trajectories  of  the  accelerometers  were 
computed  and  plotted. 

The  zero  time  for  the  trajectory  computation  was  taken  as  the 
earliest  positive  rise  on  an  accelerometer  waveform.  Apparently  in  the 
pin  withdrawal  process,  a  short  duration  negative  acceleration  can  be 
produced.  This  initial  disturbance  on  an  accelerometer  record  was  not 
used  in  the  trajectory  computations. 

B .  Tests  and  Results 

A  number  of  tests  were  conducted  in  the  test  facility  using  both  the 
Endevco  and  the  Schaevitz  accelerometers.  Table  5  lists  those  reported 
here . 

Tests  28,  29,  30,  43,  and  45  were  run  with  springs  that  produced  peak 
accelerations  of  about  4-4.5  g's.  Tests  31,  33,  35,  46,  and  48  were  run 
with  springs  that  produced  peak  accelerations  of  about  7-9  g's. 

The  results  for  tests  29,  33,  45,  and  46  will  be  presented  in  some 
detail.  Tests  29  and  33  are  low  and  high  g  tests  for  the  sensor  arrav 
with  Endevco  accelerometers,  and  45  and  46  are  low  and  high  g  tests  with 
the  Schaevitz  accelerometers. 
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TABLE  5.  TESTS  PERFORMED  WITH  INERTIAL  SENSOR  ARRAYS  IN  THE  TEST  FACILITY 


Test 

Number 


Long 

Spring 

Force 

(N) 


Total 

Spring 

Force 

(N) 


XZero* 

(m) 


Two  Endevco  accelerometers,  one  Humphrey  rate  sesnor 


Two  Schaevitz  accelerometers,  one  Humphrey  rate  sensor 


Radius 

U) 


28 

214 

712 

0.963 

1.391 

29 

214 

712 

0.964 

1.391 

30 

214 

712 

0.965 

1.391 

31 

400 

1334 

0.962 

1.391 

33 

400 

1334 

0.969 

1.391 

35 

400 

1334 

0.963 

1.391 

43 

222 

712 

0.982 

1.386 

45 

285 

801 

0.983 

1.386 

46 

400 

1334 

0.983 

1.386 

48 

400 

1472 

0.980 

1.386 

*XZero  is  the  horizontal  distance  from  the  axis  to  the  vertical  line 
through  the  sensitive  masses  of  the  accelerometers. 


This  test  was  done  using  the  Endevco  accelerometers.  The  total 
force  in  the  cable  at  initiation  of  the  test  was  712  Newtons,  and  a  peak 
acceleration  of  about  4.5  g's  was  produced. 

Figure  19  shows  the  unfiltered  vertical  acceleration  record,  and 
Figure  20  presents  the  filtered  record.  The  large  negative  deflection  at 
about  1.6  seconds  corresponds  to  rotation  of  the  quarter  wheel  90  degrees, 
where  it  impacts  a  plastic  foam  stop  and  bounces.  Figures  21  and  22  show 
the  unfiltered  and  filtered  records  for  the  initially  horizontal 
accelerometer . 

The  records  for  both  accelerometers  show  an  initial  negative 
acceleration.  Lt  was  assumed  that  this  was  produced  during  the  process  of 
releasing  the  restraint  cable  by  yanking  the  pin  out  of  the  holder.  See 
Figure  13.  This  negative  pulse  was  omitted  from  the  waveform  used  in  the 
data  reduction  program.  Zero  time  was  taken  as  the  time  for  the  first 
positive  rise  of  an  accelerometer  waveform.  This  time  was  used  as  zero 
time  for  the  other  waveforms  as  well. 

Figures  23  and  24  show  the  unfiltered  and  filtered  rate  sensor 
records,  respectively.  Because  the  rate  sensor  has  a  very  low  natural 
frequency,  it  has  a  time  lag  in  regions  of  rapidly  changing  angular 
velocities.  As  a  result,  the  waveform  initially  has  negative  values. 

Figures  25  and  26  show  the  unfiltered  and  filtered  records  for 
the  angle  of  rotation  as  derived  from  the  potentiometer  on  the  axis  of  the 
quarter  wheel.  The  bounce  of  the  wheel  is  evident. 

Figures  27  and  28  show  the  unfiltered  and  filtered  records  from 
the  load  cell.  Initially  the  peak  load  is  on  the  cable.  When  the 
restraint  cable  pin  is  pulled,  the  strong,  short  spring  unloads  rapidly. 
The  weaker,  long  spring  unloads  at  a  much  slower  rate. 

Figures  29  and  30  show  the  unfiltered  and  filtered  positive 
phases  of  acceleration  from  the  initially  vertical  accelerometer.  Figures 
31  and  32  show  the  unfiltered  and  filtered  records  for  the  initially 
horizontal  accelerometers  for  the  same  time  period.  Figures  33  and  34 
show  the  unfiltered  and  filtered  load  cell  waveforms.  As  expected,  the 
shapes  of  the  acceleration  records  are  similar  to  that  of  the  load  cell. 

Figure  35  shows  the  trajectory  of  the  accelerometers.  The  solid 
reference  curve  is  that  determined  from  the  radius  to  the  accelerometers 
ard  the  angle  from  the  potentiometer  record.  The  dashed  curve  is  that 
derived  from  the  accelerometer  and  rate  sensor  records  using  the  programs 
described  in  Section  VI  and  listed  in  Appendix  A.  Both  curves  terminate 
at  the  same  time,  which  is  the  last  time  that  the  data  are  uninfluenced  by 
deceleration  of  the  quarter  wheel  as  it  impacts  the  plastic  foam  stop. 
Figure  36  shows  the  horizontal  X  error  versus  time.  For  most  of  the  time, 
the  error  is  less  than  0.03  metres  (1.2  inches).  At  the  end  point  the 
error  is  -0.086  metres  (-3.4  inches).  Figure  37  shows  the  vertical  V 
error  versus  time.  It  remains  less  than  0.03  metres  (1.2  inches). 
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e  24.  Filtered  angular  veloci 
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Figure  31. 
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Unfiltered  positive  phase  of  acceleration  versus  time  from  the 
initially  horizontal  accelerometer  on  Test  29. 


Figure  12.  Filtered  positive  phase  of  acceleration  versus  time  from  the 
initially  horizontal  accelerometer  on  Test  29. 
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The  dashed  line  in  Figure  38  is  the  angle  of  rotation  as 
determined  by  integrating  the  filtered  rate  sensor  record.  The  solid 
curve  is  the  angle  derived  from  the  potentiometer  record.  The  agreement 
is  excellent.  Figure  39  shows  the  angular  error  versus  time.  It  remains 
less  than  one  degree  for  most  of  the  time,  and  is  always  less  than  two 
degrees  in  magnitude. 

2.  Test  33 

This  test  was  done  using  the  Endevco  accelerometers.  The  total 
force  in  the  cable  at  initiation  of  the  test  was  1334  Newtons,  and  a  peak 
acceleration  of  about  8  g's  was  produced. 

Figures  40  through  51  show  unfiltered  and  filtered  records  for 
the  accelerometers,  rate  sensor,  and  potentiometer.  The  load  cell  record 
was  bad,  and  it  is  not  included. 

Figure  52  shows  the  trajectory  of  the  accelerometers.  The  solid 
curve  is  that  determined  from  the  radius  and  angle  from  the  potentiometer 
record.  The  dashed  curve  is  that  derived  from  the  accelerometer  and  rate 
sensor  records  using  the  computer  program  listed  in  Appendix  A. 

Figure  53  shows  the  horizontal  X  error  versus  time.  Figure  54 
shows  the  vertical  Y  error  versus  time.  Both  are  0.03  metres  or  less. 
The  agreement  Is  good. 

Figure  55  compares  the  angle  of  rotation  as  determined  by 
integrating  the  rate  sensor  record  with  that  derived  from  the 
potentiometer  record.  Figure  56  shows  the  error  in  angle  versus  time. 
For  most  of  the  time,  the  error  is  less  than  2  degrees  in  magnitude.  It 
reaches  a  maximum  of  about  -3.4  degrees  at  impact. 

3.  Test  45 

This  test  was  done  using  the  Schaevitz  accelerometers.  The 
total  force  in  the  cable  at  initiation  of  the  test  was  801  Newtons,  and  a 
peak  acceleration  of  about  4.4  g's  was  produced. 

Figures  57  through  72  show  unfiltered  and  filtered  records  for 
the  accelerometers,  rate  sensor,  potentiometer,  and  load  cell.  The 
magnitude  of  the  oscillations  appearing  on  the  Schaevitz  accelerometer 
records  are  considerably  less  than  those  on  the  Endevco  accelerometer 

records . 

Figure  73  shows  the  trajectory  of  the  accelerometers.  The  solid 
curve  is  that  determined  from  the  radius  and  angle  from  the  potentiometer 
record.  The  dashed  curve  is  that  derived  from  the  accelerometer  and  rate 
sensor  records  using  the  computer  program  listed  in  Appendix  A.  Figures 
74  and  75  show  the  horizontal  and  vertical  displacement  errors.  At  impact 
they  are  -0.051  and  0.108  metres,  respectively. 

Figure  76  compares  the  angle  of  rotation  as  determined  bv 
integrating  the  rate  sensor  record  with  that  derived  from  the 
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Error  in  the  rotation  angle  obtained  by  integration  of  the  rate  sensor  record  for 
Test  33. 
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Figure  65.  Unfiltered  load  cell  record  from  Test  45 
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Figure  66.  Filtered  load  cell  record  from  Test  45. 


Figure  67.  Unf 
ver 


2/2 


TESTING  OF  AN  INERTIAL  REFERENCE  SVSTEH  CONCERT  TO 
MEASURE  8LATT-INDUCE0  (U)  TECH  REPS  INC  ABERDEEN  HD 
ABERDEEN  RESEARCH  CENTER  N  H  ETHRIDGE  ET  AL  NOV  87 
ARC-87-184  BRL-CR-591  F/G  14/2  NL 


£[tTi] 


720 


nfiltirfd 


560 

TIME,  ms 


Unfiltered  positive  phase  of  the  load  cell  record  from 
Test  45. 


Filtered  positive  phase  of  the  load  cell  record  from  Test  45. 


Figure  73.  Comparison  of  the  sensor  array  trajectory  for  Test  45  derived  from 
the  Schaevitz  accelerometers  and  the  Humphrey  rate  sensor  with  the 
reference  trajectory  derived  from  the  potentiometer  record. 


potentiometer  record.  Figure  77  shows  the  error  in  angle  versus  time. 
For  most  of  the  time,  the  error  is  less  than  2  degrees  in  magnitude.  It 

is  a  maximum  of  -2.60  degrees  at  impact. 

4 .  Test  46 

This  test  was  performed  using  the  Schaevitz  accelerometers.  The 
total  force  in  the  cable  at  initiation  of  the  test  was  1334  Newtons,  and  a 
peak,  acceleration  of  about  7.4  g's  was  produced. 

Figures  78  through  93  show  unfiltered  and  filtered  records  for 
the  accelerometers,  rate  sensor,  potentiometer,  and  load  cell.  The 
oscillations  on  the  accelerometer  waveforms  are  lower  in  magnitude  than 
those  that  occur  on  the  Endevco  accelerometer  records. 

Figure  94  shows  the  trajectory  of  the  accelerometers.  The 
dashed  curve  is  that  derived  from  the  accelerometer  and  rate  sensor 
records  using  the  computer  program  of  Appendix  A.  The  solid  reference 
curve  is  that  derived  using  the  radius  and  the  potentiometer  record.  Here 
trie  disagreement  is  large.  Figures  95  and  96  show  the  horizontal  and 
vertical  displacement  errors.  The  horizontal  displacement  error  is  small 
and  acceptable.  The  vertical  displacement  error  increases  in  time  to  an 
uLtimate  value  of  0.16  metres. 

Figure  97  compares  the  angle  of  rotation  as  determined  by 
integrating  the  rate  sensor  record  with  that  derived  from  the 
potentiometer  record.  Figure  98  shows  the  error  in  angle  versus  time. 
For  most  of  the  time  the  error  is  about  one  degree  in  magnitude.  It 
increases  to  -2.4  degrees  at  impact. 

5.  Tests  28,  30,  31,  and  35 

These  tests  were  performed  using  the  Endevco  accelerometers. 

For  Tests  28  and  30,  the  total  force  in  the  cable  at  initiation 
was  712  Newtons.  For  Tests  31  and  35,  the  force  was  1334  Newtons. 

The  peak  accelerations  were  about  4  g's  for  Tests  28  and  30,  and 
about  8.5  g's  for  Tests  31  and  35. 

Figures  99  through  118  show  the  comparisons  of  the  trajectories 
derived  from  the  inertial  sensors  with  that  derived  from  the  potentiometer 
1  ji  >  rd  and  the  radius,  the  horizontal  and  vertical  displacement  errors, 
t , i e  comparison  of  the  rotation  angle  from  the  integral  of  the  rate  sensor 
•e.  >rd  with  that  from  the  potentiometer,  and  the  error  in  angle  versus 
ir..e . 

6 .  Tests  43  and  48 

These  tests  were  performed  using  the  Schaevitz  accelerometers. 

For  Test  43,  the  total  force  in  the  cable  at  initiation  was  7!.' 
..  wr.'-ns.  The  peak  acceleration  was  about  3.5  g's. 
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Error  in  the  rotation  angle  obtained  by  integration  of  the  rate  sensor  record  for 
Test  45. 
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Figure  78.  Unfiltered  acceleration  versus  time  from  the  initially 
vertical  Schaevitz  accelerometer  on  Test  46. 
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Figure  79.  Filtered  acceleration  versus  time  from  the  initially  vertical 
Schaevitz  accelerometer  on  Test  46. 
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Figure  94.  Comparison  of  the  sensor  array  trajectory  for  Test  46  derived  from  the  Schaevitz 

accelerometers  and  the  Humphrey  rate  sensor  with  the  reference  trajectory  derived  from 
the  potentiometer  record. 


rison  of  the  rotation  angle  obtained  by  integration  of 
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Figure  99.  Comparison  of  the  sensor  array  trajectory  for  Test  28  derived  from  the  Endevco 

accelerometers  and  the  Humphrey  rate  sensor  with  the  reference  trajectory  derived  from 
the  potentiometer  record. 
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.  Horizontal  displacement  error  versus  time  of  the  sensor  array 
for  Test  28. 


Figure  101.  Vertical  displacement  error  versus  time  of  the  sensor  array 
for  Test  28. 


ANGULAR  ERROR 


obtained  by  integration 


Comparison  of  the  sensor  array  trajectory  for  Test  30  derived  from  the  Endevco 
accelerometers  and  the  Humphrey  rate  sensor  with  the  reference  trajectory  derived 
the  potentiometer  record. 


Comparison  of  the  rotation  angle  obtained  by  integration  of  the  rate  sensor  record 
that  derived  from  the  potentiometer  record  for  Test  30. 
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Figure  109.  Comparison  of  the  sensor  array  trajectory  for  Test  31  derived  from  the  Endevco 

accelerometers  and  the  Humphrey  rate  sensor  with  the  reference  trajectory  derived  from 
the  potentiometer  record. 
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jure  112.  Comparison  of  the  rotation  angle  obtained  by  integration  of  the  rate  sensor  record  with 
that  derived  from  the  potentiometer  record  for  Test  31. 
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Figure  1 11.  Comparison  of  the  rotation  angle  obtained  by  integration  of  the  rate  sensor  record  with 
that  derived  from  the  potentiometer  record  for  Test  43. 
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Comparison  of  the  rotation  angle  obtained  by  integration  of 
that  derived  from  the  potentiometer  record  for  Test  48, 


computed  angle  versus  time  for  a  truck-shelter  system  struck  by  blast 
i  Test  28  in  the  test  facility. 
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:!■  i '  lists  the  errors  in  the  horizontal  and  vertical  displacements 
•  ■  end  . ’ t  the  trajectories  in  the  test  facility.  Also  listed  is  the 
ui  the  rotation  angle  derived  from  the  integral  of  the  rate  sensor 


.he  error  in  the  rotation  angle  at  the  end  of  the  trajectory  ranges 
- 1  .  h  to  -  '.hi  degrees.  The  error  could  be  reduced  by  correcting  for 
the  I  a  r  espouse  time  that  occurs  where  the  angular  velocity  is 

.  i ::  i  i>  c  r an  i  d  1  v  . 

.ere  -.eems  to  be  no  distinction  in  errors  between  high  and  low  g 
to  -t->.  ’  e  schaevitz  accelerometers  produced  X  errors  that  are  about  the 

•  .e  the  Fndevco  X  errors.  However,  the  Schaevitz  Y  errors  are  larger 

t  h  :r.  t  he  hade  \'co  V  errors.  The  outputs  from  both  the  vertical  and 
horizontal  accelerometers  are  used  to  obtain  X  or  to  obtain  Y.  At  the  end 
o:  the  contract  effort,  no  explanation  had  been  found  for  these  larger 
errors.  The  Schaevitz  accelerometer  has  a  much  lower  natural  frequency 
tnan  the  Fndevco  accelerometer,  but  it  is  rated  for  higher  accuracy. 

i  lie  errors  from  the  tests  using  the  Endevco  accelerometers  seem  small 
enough  to  warrant  further  pursuit  of  an  inertial  reference  system  for  use 
on  vehicles.  The  records  in  the  tests  reported  here  were  obtained  using 
: -the-shel f  sensors,  and  were  processed  using  a  simple  time-marching 
integration  program.  The  accuracy  of  performance  of  a  sensor  array  should 
improve  if  filters  are  matched  to  gage  and  input  waveform  characteristics, 
and  if  the  distortions  in  the  data  produced  by  gage-response 

characteristics  are  removed  or  reduced. 

It  seems  unlikely  that  a  sensor  array  can  track  a  vehicle  struck  by 
blast  beyond  first  impact  with  the  ground,  whether  the  impact  is  due  to 
overturning  or  to  the  return  to  ground  after  the  vehicle  has  been 

airborne . 


The  sensor  array  will  require  packaging  so  that  the  arrival  of  the 
blast  wave  at  the  vehicle  does  not  disable  the  accelerometers  or  the  rate 
sensors  for  more  than  a  few  milliseconds,  if  at  all.  The  test  facility  at 
present  does  not  provide  an  impact  at  the  start  of  rotation  corresponding 
to  that  experienced  by  a  vehicle  struck  by  blast. 

VIII.  REQUIREMENTS  FOR  AN  INERTIAL  REFERENCE  SYSTEM 

."he  results  of  the  test  program  indicate  that  the  array  of  off-the- 
shelf  sensors  can  track  movement  in  a  plane  within  acceptable  error 
bounds.  Figure  131  shows  a  block  diagram  of  a  system  for  use  on  a  vehicle 
with  six  degrees-of-f reedom  of  movement. 

i ne  array  of  accelerometers  and  rate  sensors  is  packaged  separately 
r  i:  i  :ri’  tlie  size  of  the  array  and  facilitate  its  mounting  near  the 

•  •  ?  cravltv  of  the  vehicle. 


TABLE 


6.  ERRORS  IN  HORIZONTAL  X  AND  VERTICAL  Y  DISPLACEMENTS  AND  ROTATION 
ANGLE  9  AT  END  OF  TRAJECTORIES  IN  THE  TEST  FACILITY 


Test  Number 

X  Error 

Y  Error 

6  Error 

(m) 

(m) 

(deg) 

Two  Endevco 

accelerometers , 

one  Humphrey  rate 

sensor 

28 

-0.018 

0.037 

-3.14 

29 

-0.086 

0.025 

-1.81 

30 

-0. 145 

0.085 

-2.20 

31 

0.036 

0.068 

-1.20 

33 

0.023 

0.031 

-3.37 

35 

-0.135 

0.067 

-3.51 

Two  Schaevitz  accelerometers,  one  Humphrey  rate  sensor 


Block  diagram  of  inertial  reference  system 
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Hie  system  requires  an  arming  signal  prior  to  the  test  so  that  the 
equipment  tan  warm-up  and  stabilize.  This  signal  can  he  provided  either 
by  internal  programming  or  hard  wire  from  a  sequence  timer  or  manual 
s  -  itch. 


A  means  of  triggering  the  data  storage  mode  at  the  desired  time  is 
required.  The  pressure  sensor  provides  this  signal  when  the  blast  wave 
arrives  at  its  location.  The  sensor  is  to  be  placed  on  or  near  the  outer 
surface  of  the  vehicle  facing  the  blast  so  that  the  data  recording  will 
commence  before  any  sensor  response  takes  place. 


Once  stored  in  memory,  the 
retrieval  process  takes  place, 
sequence  to  the  computer  where  it 
disk. 


data  must  remain  unchanged  until  the 
At  this  time  the  data  is  transferred  in 
is  stored  permanently  on  a  magnetic 


An  examination  of  vehicle  overturning  calculations  that  were  made 
using  the  BLOM  code  showed  that  a  2  1/2  ton  truck-shelter  system,  an  Ml  13 
personnel  carrier,  and  the  Infantry  Fighting  Vehicle  can  require  from  2.0 
to  3.1  seconds  to  rotate  to  an  angle  at  which  they  almost  overturn.  These 
times  applied  over  the  yield  range  for  a  nuclear  explosion  from  1  KT  to 
1000  KT.  The  time  for  fall  back  to  the  original  position  or  to  overturn 
90  degrees  will  be  approximately  the  same.  Therefore,  to  track  the  full 
rotational  movement  of  a  vehicle  under  conditions  where  it  is  exposed  near 
tiie  threshold  for  overturning,  a  recording  time  of  at  least  seven  seconds 
is  required. 


The  data  from  the  gages  are  filtered  before  they  are  digitized.  The 
filters  and  sampling  rates  are  selected  considering  the  frequency  content 
of  the  sensor  signals.  The  records  from  the  Endevco  and  Schaevitz 
accelerometers  used  in  the  tests  described  in  this  report  still  retained 
significant  oscillations  at  a  frequency  of  about  250  Hz  after  passing 
through  a  100  Hz  low-pass  filter.  For  reducing  the  error  due  to  sampling 
rate  to  about  one  percent,  the  sampling  rate  should  be  5000  Hz  or  more. 

The  minimum  storage  capacity  for  each  channel,  considering  seven 
seconds  of  recording  at  a  rate  of  5000  Hz,  is  35000  data  points.  With  a 
storage  capacity  of  64000  12  bit  words,  eight  seconds  of  recording  could 
he  obtained  at  a  sampling  rate  of  8000  Hz.  Such  a  sampling  rate  is  fast 
enough  to  ensure  that  no  information  will  be  lost  due  to  the  sampling 
:>  roc  ess.  The  sampling  times  are  synchronized  for  all  channels. 

The  system  is  computer  driven.  The  computer  is  used  to  give 
.!  i  :ve l ions  to  tne  system,  insert  calibration  data,  monitor  system  status, 
v ■<!  retrieve  data. 

For  use  on  a  large-scale  blast  o  .emulation  in  the  field,  the  system 
1  ■  be  capable  of  operating  unattended  for  a  period  of  7  2  hours,  where 
a-::  t  :  "it  ion  'ho  a  ->  ingle  event  may  occur  within  the  first  48  hours,  and 
•  t  •  r.  c  not  h  •  ret-ieved  f  ,>r  another  24  hours. 


Kv 


i'  V  lists  characteristics  of  the  system.  There  are  several 
verm  <•  re  i.ii  Lv  available  systems  that  may  meet  the  hardware  requirements 
a: ter  some  mod  if ieations.  However,  the  software  must  be  developed. 

XI.  CONCLUSION’S 

Alter  a  review  ot  the  literature,  the  conclusion  was  made  that  the 
best  combination  of  sensors  for  an  inertial  reference  system  for  tracking 
tin-  movement  of  a  vehicle  struck  by  blast  is  three  linear  accelerometers 

ill  t  rce  angular  rate  sensors.  This  conclusion  was  made  assuming  that 

.ages  with,  satisfactory  response  characteristics  can  be  found. 

before  proceeding  with  building  a  complete  six  degrees-of-f reedom 
inertial  reference  system,  it  is  highly  desirable  to  verify  the  tracking 
capability  of  sensor  arrays  for  a  simpler  system. 

It  was  found  possible  to  construct  a  test  facility  that  would 
simulate  tile  overturning  of  a  vehicle  struck  by  blast.  The  required 
sensor  array  tor  tracking  rotation  in  a  plane  is  two  linear  accelerometers 

and  one  angular  rate  sensor. 

Two  arrays  were  assembled  and  tested  on  the  test  facility.  One  used 
Endevco  accelerometers  and  the  other  used  Schaevitz  servo  accelerometers. 
After  testing  the  sensor  arrays  in  the  test  facility,  it  was  concluded 

that  tiie  sensor  array  with  off-the-shelf  Endevco  accelerometers  and  a 

Humphrey  rate  sensor  tracked  displacement  and  rotation  in  a  plane  within 
acceptable  error  bounds. 

Other  than  filtering  with  a  100  Hz  low-pass  filter,  no  preprocessing 
or  modification  of  the  gage  records  was  performed.  The  accuracy  of 
performance  of  the  sensor  array  can  be  improved  if  filters  are  matched  to 
gage  and  input  waveform  characteristics,  and  if  the  distortions  in  the 
data  produced  by  gage  response  characteristics  are  removed  or  reduced. 

The  results  of  the  tests  show  that  the  development  of  a  complete 
inertial  reference  system  is  worth  pursuing.  The  general  requirements  for 
such  a  system  were  determined. 

Because  of  the  large  signals  generated  on  the  gage  records  by  impact 
oi  the  venicle  with  the  ground,  an  inertial  reference  system  cannot  be 
expected  to  track  the  vehicle  after  the  first  impact  with  the  ground. 

me  sensor  array  with  the  Schaevitz  accelerometers  did  not  track 
d  i  :-p :  ;k  emenr  satisfactorily.  At  the  time  of  the  end  of  the  contract,  no 
•  :  an.it  ion  for  the  large  error  had  been  found.  The  Schaevitz 

acre  1 e  rone  ter  has  a  much  lower  natural  frequency  than  the  Endevco 
.n  ,  o  1  »-r  c  r  ,  but  is  rated  for  higher  accuracy. 

•  ■  ■  ir : 1  it v  was  simple  to  use  and  performed  well.  The  means  ot 

V:  : he  restraint  cable  needs  to  be  improved  to  minimize  the 
1 :  ;  e  ■  the  waveforms. 


TABLE  7.  CHARACTERISTICS  OF  INERTIAL  REFERENCE  SYSTEM 

FOR  TRACKING  THE  BLAST  OVERTURNING  OF  VEHICLES 


Physical  Construction 

*  Packaged  as  two  components  connected  by  electrical  cable.  One 
component  is  the  sensor  array,  the  other  is  the  recording  and 
power  package. 

*  Sensors  mounted  for  mechanical  filtering  of  vibration  induced  by 
blast 

*  Both  components  of  rugged,  shock-  and  vibration-resistant 
construction 


Gages 


*  Three  linear  accelerometers  (Endevco) 

*  Three  rate  sensors  (Humphrey,  Inc.) 


Memories 


*  Non-volatile 


*  Minimum  64K  12  bit  words  per  channel 


A/D  Converters 


*  12  bit 


D 


<vc 


*  4000  and  8000  Hz  sampling  rates 

*  Synchronized 
Operating  Modes 

*  Activated  by  manual  switch,  sequence  timer,  or  internal  clock- 
controlled  program 

*  Armed  by  manual  switch,  sequence  timer,  or  internal  clock- 

controlled  program 


Triggered  for  data  storage  by  pressure  sensor  signal,  manual 
switch,  sequence  timer,  or  internal  clock-controlled  program 


Control  Unit 


*  Controls  operation  of  system  components 

*  Receives  external  signals  and  instructions 

*  Has  capability  for  internal  clock-driven  program 
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TABLE  7.  CHARACTERISTICS  OF  INERTIAL  REFERENCE  SYSTEM  FOR 

TRACKING  BLAST  OVERTURNING  OF  VEHICLES  (CONTINUED) 


Cal ibra  t ion 


Control  unit  directed 


Use  of  shunt  resistor  for  accelerometer  channels 


Use  of  voltage  substitution  for  rate  sensor  channels 


Power  Supply 


Self-contained,  rechargeable  batteries 


Computer 


Adequate  to  support  system  for  72  hours  for  a  single  event  test, 
with  the  single  event  occurring  in  the  first  48  hours 

Adequate  to  maintain  memory  for  months 


Delivers  instructions  to  control  unit 


*  Interrogates  system  to  monitor  status 

*  Reads  data  from  system  and  stores  data  on  disks 

*  Processes  records  to  remove  noise  and  correct  for  sensor 
d  istortions 


Computes  movement  of  vehicle 


Sof  tware 


Available  to  perform  filtering,  noise  removal,  spectral  analysis 
and  waveform  reconstitution,  and  corrections  for  sensor- 
introduced  distortions 


*  Program  to  utilize  corrected  records  to  derive  displacement  and 
orientation  versus  time  of  the  vehicle 
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APPENDIX  A 

DATA  REDUCTION  PROGRAM 


The  following  program  was  written  and  compiled 
on  the  IBM  PC  AT  using  version  1.10  of  the  Ryan- 
McFarland  RM/ FORTRAN  language  compiler. 


£*.#*#♦♦«>«►♦'**»**•****************♦**♦******♦******■*****♦*  ****************' 

c 

c 

C  TEST  FIXTURE  DATA  REDUCTION 

C  ABERDEEN  RESEARCH  CENTER 

C  NOEL  H.  ETHRIDGE 

C 
C 

£***«.  ****«*»*****»***«►*  ********«**»***  •***«»»«•**»*»  »****•****•*•****••• 

c 

C  PROGRAM  DESCRIPTION 

C 

C 

C  THIS  PROGRAM  IS  INTENDED  TO  COMPUTE  THE  UERTICAL  AND 

C  HORIZONTAL  MOUEMENT  OF  A  SENSOR  ARRAY  ROTATING  IN  A  UERTICAL 

C  PLANE  ABOUT  A  FIXED  AXIS.  THE  SENSOR  ARRAY  IS  MOUNTED  ON  A 

C  FIXTURE  THAT  IS  ESSENTIALLY  A  WHEEL.  THE  WHEEL  IS  LOADED  BY 

C  A  CABLE  ATTACHED  TO  ITS  PERIPHERY  AND  HELD  IN  PLACE  BY  ANOTHER 

C  CABLE.  AT  START  THE  RESTRAINT  CABLE  IS  SUDDENLY  RELEASED  AND 

C  THE  WHEEL  ROTATES  THROUGH  APPROXIMATELY  SO  DEGREES.  THE  SENSOR 

C  ARRAY  MOUES  THROUGH  A  KNOWN  ARC.  THE  OBJECT  OF  THE  TEST  AND 

C  THIS  COMPUTATION  IS  TO  COMPARE  THE  TRAJECTORY  OF  THE  SENSOR 

C  ARRAY  AS  DETERMINED  FROM  THE  SENSOR  RECORDS  WITH  THAT  DETERMINED 

C  BY  OTHER  MEANS,  SUCH  AS  A  COMPUTATION  USING  THE  MEASURED  ANGLE 

C  OF  ROTATION  AND  THE  RADIUS  FROM  THE  ROTATION  AXIS  TO  THE  SENSOR 

C  ARRAY . 

C 

C  THE  ARRAY  CONSISTS  OF  TWO  LINEAR  ACCELEROMETERS  AND  ONE 

C  ANGULAR  RATE  SENSOR.  ONE  ACCELEROMETER  IS  INITIALLY  UERTICAL 

C  AND  THE  OTHER  HORIZONTAL.  OTHER  SENSORS  USED  ARE  A  LOAD  CELL  TO 

C  MEASURE  THE  FORCE  IN  THE  CABLE  AND  A  PRECISION  POTENTIOMETER 

C  ATTACHED  TO  THE  ROTATION  AXIS  TO  MEASURE  ROTATION  ANGLE. 

C 

C  THE  ANGULAR  UELOCITY  FRDM  THE  RATE  SENSOR  IS  INTEGRATED  TO 

C  OBTAIN  THE  ANGLE  OF  ROTATION.  THIS  ANGLE  IS  USED  TO  CORRECT  THE 

C  READINGS  OF  THE  ACCELEROMETERS  FOR  THE  GRAUITY  FIELD,  AND  TO 

C  CALCULATE  THE  HORIZONTAL  X  AND  UERTICAL  Y  COMPONENTS  OF  ACCELER- 

C  AT  I  ON  FROM  THE  CORRECTED  ACCELEROMETER  RECORDS.  THE  X  AND  Y 

C  COMPONENTS  OF  ACCELERATION  ARE  DOUBLE  INTEGRATED  TO  OBTAIN 

C  X  AND  Y  UELOCITIES  AND  DISPLACEMENTS. 

C 

C  AN  ALTERNATIUE  COMPUTATION  IS  PERFORMED  IN  WHICH  THE  ANGLE 

C  TO  CORRECT  THE  GRAUITY  FIELD  AND  TO  DETERMINE  THE  X  AND  Y  COMPO- 

C  NENTS  IS  OBTAINED  FROM  THE  POTENTIOMETER  RECORD.  THE  COMPUTED 

C  MOUEMENT  DERIUED  USING  THE  ROTATION  ANGLE  MEASURED  BY  THE 

C  POTENT  I OMENTER  AND  THE  RADIUS  TO  THE  SENSOR  IS  PRESENTED  AS  A 

C  REFERENCE  FOR  EUALUATION  OF  THE  RESULTS  DERIUED  FROM  THE  SENSOR 

C  ARRAY.  LAST  RECORD  IS  SSS3  S99S  9999  9999  9999. 

C 


c 


DEFINITIONS  OF  PARAMETERS 


C 

C 


]  r! 


c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 
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c 

c 

c 
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C 
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AT1TLE  «  IDENTIFICATION  OF  TEST,  DATE 
GAGES  -  GAGE  DATA 


SPRNGS 

TSTEP 

NSTOP 

RADIUS 

X2ER0 

A2ER0 

A2ER0D 

G 

TIDE 

ANGUEL 

ACH 

ACU 

ANGPD 

FORCE 

ANGRSD 

ANGRS 

ACHCG 

ACUCG 

ACXHU 


ACYHU 


UELACX 

UELACY 

XACRS 

YACRS 


XRPOT 


YRPOT 

ANGPOT 

ACMCGP 


AOUCGP 
AC  X  HUP 


SPRING  DATA 

TIME  STEP  -  SAMPL I NG  INTERUAL  CS) 

NUMBER  OF  TIME  STEP  TO  STOP  COMPUTATIONS 

RADIUS  FROM  AXIS  OF  ROTATION  TO  AXES  OF  ACCELEROMETERS  (M) 
HOP  1 20NTAL  DISTANCE  FROM  AXIS  TO  ACCELEROMETER  AXES  CM) 
INITIAL  ANGLE  OF  RADIUS  WITH  H0RI20NTAL  C RAD IANS) 

A2ER0  EXPRESSED  IN  DEGREES  CDEG) 

ACCELERATION  DUE  TO  GRAUITY  -  3 . BOO  M/S**E  AT  ABERDEEN 
TIME  FROM  START  OF  LOADING  CS) 

ANGULAR  UELOCITY  FROM  RATE  SENSOR  CDEG/S) 

ACCELERATION  MEASURED  BY  ACCELEROMETER  WITH  AXIS 
INITIALLY  H0RI20NTAL  CG'S) 

ACCELERATION  MEASURED  BY  ACCELEROMETER  WITH  AXIS 
INITIALLY  UERTICAL  CG’S) 

ANGLE  OF  ROTATION  FROM  POTENTIOMETER  CDEG) 

FORCE  IN  LOADING  CABLE  MEASURED  BY  LOAD  CELL  CNEWTONS) 
ANGLE  ROTATED  FROM  START  DERIUED  BY  INTEGRATING  RATE 
SENSOR  RECORD  CDEG) 

ANGRSD  EXPRESSED  IN  RADIANS  C RAD IANS) 

ACCELERATION  MEASURED  BY  HORIZONTAL  ACCELEROMETER 
CORRECTED  FOR  GRAUITATIONAL  FIELD  CG’S) 

ACCELERATION  MEASURED  BY  UERTICAL  ACCELEROMETER 
CORRECTED  FOR  GRAUITATIONAL  FIELD  CG’S) 

X  COMPONENT  OF  ACCELERATION  OBTAINED  BY  COMBINING  X 
COMPONENTS  FROM  BOTH  CORRECTED  ACCELERATION  RECORDS 
CM/S**2) 

Y  COMPONENT  OF  ACCELERATION  OBTAINED  BY  COMBINING  Y 
COMPONENTS  FROM  BOTH  CORRECTED  ACCELERATION  RECORDS 

CM/S*«S) 

FIRST  INTEGRAL  OF  ACXHU  CM/S) 

FIRST  INTEGRAL  OF  ACYHU  CM/S) 

X  DISPLACEMENT  OF  SENSOR  ARRAY  DERIUED  USING  ANGLE  FROM 
RATE  SENSOR,  SECOND  INTEGRAL  OF  ACXHU  CM) 

Y  DISPLACEMENT  OF  SENSOR  ARRAY  DERIUED  USING  ANGLE  FROM 
RATE  SENSOR,  SECOND  INTEGRAL  OF  ACYHU  CM) 

X  DISPLACEMENT  OF  SENSOR  ARRAY  DERIUED  USING  ANGLE  FROM 
POTENTIOMETER  AND  RADIUS  TO  SENSOR  ARRAY  CM) 

Y  DISPLACEMENT  OF  SENSOR  ARRAY  DERIUED  USING  ANGLE  FROM 
POTENTIOMETER  AND  RADIUS  TO  SENSOR  ARRAY  CM) 

ANGPD  EXPRESSED  IN  RADIANS  C RAD IANS) 

HORIZONTAL  ACCELEROMETER  READING  ACH  CORRECTED  FOR 
GRAUITATIONAL  FIELD  USING  ANGLE  FROM  POTENTIOMETER  CG’S) 
UERTICAL  ACCELEROMETER  READING  ACU  CORRECTED  FOR  GRAUI¬ 
TATIONAL  FIELD  USING  ANGLE  FROM  POTENTIOMETER  CG’S) 

X  COMPONENT  OF  ACCELERATION  OBTAINED  BY  COMBINING  X 
COMPONENTS  FROM  BOTH  CORRECTED  ACCELERATION  RECORDS 
-HERE  CORRECTION  WAS  MADE  USING  ANGLE  FROM  POTENTIOMETER 
r  m  -  o  *  *  ?  i 


Y  COMPONENT  OF  ACCELERATION  OBTAINED  BY  COMBINING  Y 
COMPONENTS  FROM  BOTH  CORRECTED  ACCELERATION  RECORDS 
WHERE  CORRECTION  WAS  MADE  USING  ANGLE  FROM  POTENTIOMETER 

,  i 


A 
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C  VELAXP 
C  VELAYP 
C  XACPOT 

C  YACPOT 


C 

c 

C  MAIN  PROGRAM 

C 

c 

DIMENSION  TIMEC4000)  ANGUELC4000) ,  ACHC4000),  ACUC4000) , 

1  ANGPDC4000) ,  F0RCEC4000),  ANGRSDC2) ,  ANGRSC2) ,  ACVCGC2),  ACHCGC2) 
£  ,  ACXHVC2),  ACYHVC2),  VELACXC2) ,  VELACYC2),  XACRSC2),  YACRSC2) , 

3  XRPOT ( 4000 ) ,  YRPOT C 4000 ) ,  ANGP0TC2) ,  ACHCGPC2) ,  ACUCGPC2) , 

4  ACXHVPC2) ,  ACYHVPC2) ,  VELAXPC2),  VELAYPC2) ,  XACPQTC2) ,  YACP0TC2) 

CHARACTER*80  ATITLE 
CHARACTER*BO  GAGES 
CHARACTER*80  SPRNGS 
PI  -  3.141592654 
DGTORAD  -  PI /ISO. 

READ  IN  TITLE  OF  TEST 

READC5.il  ATITLE 
1  FORMAT  CABO) 

READ  IN  GAGE  DATA 

READ  C5.1)  GAGES 

READ  IN  SPRING  DATA 

READ  C5 , 1 )  SPRNGS 

READ  IN  RADIUS  AND  XZERO 

READ  C5.2)  RADIUS,  XZERO 
£  FORMAT  CF10.5.F11 .4) 

C 

C  READ  TIME  STEP  -  SAMPLING  INTERVAL,  AND  NUMBER  OF  TIME  STEP  FOR 

STOPPING  COMPUTATION 

C 

READ  C  5 , 3 )  TSTEP ,  NSTOP 
3  FORMAT  CF12.B, 15) 

C 

C  READ  IN  ACCELEROMETER,  RATE  SENSOR,  POTENTIOMETER,  AND  FORCE 

V  RECORDS 

DO  5  1=1 ,  NSTOP 

READ  C 5 , 4 )  ACHCI),  ACUCI),  ANGUELCI),  ANGPDCI),  FORCECI) 
t  FORMAT  C5F10.4) 

IF  CACHr I) .E0.9999)  GOTO  £0 


=  FIRST  INTEGRAL  OF  ACXHUP  CM/S) 

-  FIRST  INTEGRAL  OF  ACYHUP  CM/S) 

-  X  DISPLACEMENT  OF  SENSOR  ARRAY  DERIVED  USING  ANGLE  FROM 
POTENTIOMETER  AND  ACCELEROMETER  RECORDS  CM) 

-  Y  DISPLACEMENT  OF  SENSOR  ARRAY  DERIVED  USING  ANGLE  FROM 
POTENT  I OMENTER  AND  ACCELEROMETER  RECORDS  CM) 


c 

c 

c 

c 


c 

c 

c 

c 

100 


c 

r» 

u 

c 

c 


c 

c 

c 

c 


TIME (I)  =  FL0ATCI-1)  *  TSTEP 
FORCECI)  =  FORCE C 1 1  *  4.44B222 
CONTINUE 
NSTOP  =1-1 

DEFINE  GRAUITATIONAL  CONSTANT  Cn/S**21  AT  ABERDEEN 
G  =  3.800 

COMPUTE  INITIAL  ANGLE 

A2ER0  -  ACQSC  XZERO/RAD I  US ) 

AZEROD  =  A2ER0  *  Cl./DGTORAD) 

DEFINE  INITIAL  UALUES 

ANGRSC 1 )  =  0. 

ANGRSDC 1 )  =  0. 

UELACXC1)  =  0. 

UELACYC 1 )  =  0. 

XACRSC 1 )  -  0. 

YACRSC1)  =  0. 

UELAXPC 1 )  “  0. 

UELAYPC1)  =  0. 

XACPOT Cl)  «  0. 

YACP0TC1)  =  0. 


COMPUTE  START  UALUES  ASSUMING  ACCELEROMETERS  ARE  HORIZONTAL  AND 
UERTICAL 

ACHCGC1)  -  ACH Cl) 

ACUCGC1)  =  ACU  Cl) 

ACXHUrn  -  CACUCGC 1 )*SINCO  .  )  +  ACHCGC 1 )*COSCO . ))*G 
ACYHUC1)  =  f ACUCGC 1 )*COSCO . )  -  ACHCGC 1 )*5INC0 . ) )»G 

COMPUTE  SECOND  STEP  UALUES  AND  PERFORM  INTEGRATIONS, 

INTEGRATE  ANGULAR  UELOCITY 

ANGRSD C  2  )  -  0 . 5*C ANGUELC I +1 1  +  ANGUELC I) )*CTIMEC 1  +  1 5-TIMEC II )  «■ 
1  ANGRSDC 1) 

ANGRSC2)  =  ANGRSDC  2 1 *DGTQRAD 

CORRECT  ACCELEROMETER  RECORDS  FOR  GRAUITY  FIELD  USING  ANGLE 
OBTAINED  BY  INTEGRATING  ANGULAR  UELOCITY 

ACHCG C  2 )  -  ACHCI  +  1)  +  S INC ANGRSC 2 ) ) 

ACUCG C 8 )  =  ACU C I + 1 )  +  1.  -  COS C ANGRSC 2) ) 

COMPUTE  COMBINED  X  AND  Y  ACCELERATIONS  FROM  ACCELEROMETER 
RECORDS  USING  ANGLE  OBTAINED  BY  INTEGRATING  ANGULAR  UELOCITY 

ACXHUC  2 )  =  f ACUCG C  5 I*SINC ANGRSC 2 ) )  «■  ACHCG C 2 )*COSC ANGRSC 2 ) n*G 
ACYHU  c  2 1  =  C  ACUCG  C2)*C0SC  ANGRSC  2 ) )  -  ACHCGC2)*SINC ANGRSC2) ) )*G 


INTEGRATE  X  AND  V  ACCELERATIONS  TO  OBTAIN  X  UELDCITY  AND  DISPLAC 
r 1 F N T  AND  Y  UFLOCITY  AND  DISPLACEMENT 


UE'.L  ACX  C  E  '  -  G.‘i*«.ACXHU(2i  «■  ACXHU  ■  I  )  j  *  f  1 1  ME  C  I  + 1 1  -I  IME  <\  )  )  + 
vELACXC  x  ' 

UELACY  (  0  )  “  0 .5*CACYHUCE)  +  ACYHUt  1 )  )  *  C  T  I ME  C  I  + 1  ) -T I  MEC  I  )  )  + 


XACRS'E)  -  0 ,5*iUELACXCE)  +  UELACX C 1 ) ) * C T I  ME C l +1 ) - T IMEC I ) )  + 
1  XACRS*  1  ') 

YACRS  C  E  )  0 . S-CUELACYCE)  +  UELACYC1  )  ^CTIMECI  +1  )-TIMEC  I)  )  + 

1  YACRS C 1  i 


REPEAT  COMPUTATION  USING  ANGLE  FROM  POTENTIOMETER  RECORD 


IF  C 1  . GT . 1 5  GO  TO  105 
ANGPOTfl)  -  ANGPDC 1 )*DGTQRAD 


CORRECT  ACCELERATION  RECORDS  FOR  G  FIELD 


ACHCGPC1)  **  ACH  Cl)  +  SINC  ANGPGT  Cl  )  ) 
ACUCGPC 1)  -  ACUC1)  +  1.  -  COSC ANGPOTC 1 ) ) 


COMPUTE  COMBINED  X  AND  Y  ACCELERATIONS  FROM  BOTH  ACCELEROMETER 
RECORDS  USING  ANGLE  FROM  POTENTIOMETER 


ACXHUPC1)  -  C  ACUCGPC 1 )*S I NC ANGPOT Cl))  +  ACHCBPCl)* 
1  COS  (  ANGPOT  Cl))) *G 

ACY'HUPC  1  )  -  C  ACUCGPC  1  )*C05C ANGPOT  Cl))  -  ACHCGPC1)* 
I  S I NC  ANGPOT ( 1 ) ) ) *G 


COMPUTE  SECOND  STEP  UALUES  AND  PERFORM  INTEGRATIONS 


ANGPOT r  P '  -  ANGPOC I  +  1)*DGT0RAD 


CORRECT  ACCELEROMETER  RECORDS  FOR  GRAUITY  FIELD 


ACHCGPCE)  -  ACHCI+1)  +  SINC  ANGPOT C  E  )  ) 
ACUCGP t E )  -  ACUCI+1)  ♦  1.  -  COS C ANGPOT C P 1 > 


COMPUTE  COMBINED  X  AND  Y  ACCELERATIONS  FROM  BOTH  ACCELEROMETER 
RECORDS  USING  ANGLE  OBTAINED  FROM  POTENT  I OMETFR 


ACXHUPCE)  -  (  ACUCGP  C  E  5  *S  INC  ANGPOT  C  Em  -*•  ACHCGPCE)* 

1  COSC ANGPGT CE? ) )*G 

ACYHUPCE)  -  C  ACUCGPC  E  )*CDSC  ANGPGTC?)  )  -  ACHCGPCE)* 
1  SINC ANGPOT CE) ) )*G 


INTEGRATE  X  AND  Y  ACCELERATIONS  TO  OBTAIN  X  UELDCITY  AND 
DISPLACEMENT  AND  Y  UELDCITY  AND  DISPLACEMENT 


UELAXPi'E)  -  0  .S»C  ACXHUPCE)  +  ACXHUP  C  1  )  '  *  <  T  I  MF.C  1  ♦  1  ) -T  I  ME  C  I  )  )  + 

1  UELAXF’i'  ]  t 

UELAVPC  I?  >  -  0 . 5*  C  ACYHUPC  E  )  +  ACYHUpf  l  )  .*'7  IMF-'  I  ♦  1  liMFC  I  )  )  ♦ 
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UELAYPC  i  ) 

XftCPiJTi?)  =  0 ,5*(UELAXPC2)  +  UELAXPC  1 5  5 *C T IMEC  I  +1 5 -T IMEC  I  )  5  + 

XACPOTC 1  i 

YACP0TC25  =  C .5* (UELAYPC 2 5  +  UELAYPC  1  5  5*  C  TIMEC  I  +1  )-T IMEC  1)5  * 

YACPOT :  1 

CALCULATE  X  AND  Y  DISPLACEMENTS  FROM  RADIUS  AND  ANGLE 
IF  C !  .GT  .  1 5  GO  TO  10G 

XRPOI(l)  »  RAD 1US*  CCOSC  A2ERD 5 -COSC  C  ANGPDC 1 ) *DGTORAD) 

--  A2ERO  5  5 

YRPDT (  1  .'i  «  RADIUS* C  S I NC  C  ANGPDC  1  5*DGTORAD5  +  AZERO)  - 

SINC  A2ERO 5  5 

XRPGTCI+1)  -  R AD  I US *C CDS C AZERO)  -  COSC C ANGPDC 1+1 )*CPI / 

ISO .55  +  A2ER0 ) 5 

YRPOTCI+1)  -  RAD  I  US* CS INC  C ANGPDC 1  +  1 )*DGTORAD)  +  AZERO) 

-  SINCA2ER0) ) 

PRINT  RESULTS 

I FC I .GT . 1 5  GO  TO  101 
WRITE  16,6) 

FORMAT  C1H15 
WRITE  C  S , 7 1 
FORMAT  C/5 

WRITE  (6,05  ATITLE 
FORMAT  C  5X , ABO  5 
WRITE  (6,75 
WRITE  C6,B)  GAGES 
WRITE  (6,75 
WRITE  15,05  SPRNGS 
WRITE  C  6 , 7 5 
WRITE  C 6 , 7 5 
WRITE  C6.S5  RADIUS 

FORMAT  C10X,  SHRADIUS  -  ,  F10.4,  BH  METRES5 

WRITE  (6,105  X2ER0 

FORMAT  C 1 IX ,  0HX2ERO  -  .F10.4,  0H  METRES) 

WRITE  (6,115  A2ER0D 

FORMAT  (10X,  9HA2ER0D  -  .F10.4,  SH  DEGREES) 

WRITE  rB, 155  G 

FORMAT  C15X,  4HG  -  ,  F10.4,  BH  M/S**2) 

WRITE  (.5,6  5 
WRITE  15,71 
W R i  I E  -6,135 

FORMAT  C 1 >  BH  TIME  , 4X , 5HANGPD , 2X , 6HANGRSD , 2X , BHANGUEL , 4X , 3HACH , 

5/ ,  iHACu , 3X , 5HF0RCE , 4X , SHACHCG , 3X , 5HACUCG , 3X , 5HACXHU , 3X , 5HACYHU , 
0 w  SHXACRS  IX, 6HXACP0T , 2X , 5HXRP0T , 2X , 5HYACRS , IX , BHYACPOT , 2X , 

VaHYPPOT : 

-R IT1  -6,145 

6',:-  ’A-  .  •  ,  *H  IS;  ,  6X,5H(DEG5  ,PX,  SHC  DEG  5 , 2X  ,  7H  C  DEG/S  5 , 3X  ,  5HC  G  ’55, 

MX  SH'  G  51 , 3X , 3HCN5 , 5X , 5HCG '55 , 3X,5HCG’S) , IX , 7HCM/S*»2 , IX , 

M/f.-**P  3X.3HCM)  ,4X,3H(M5,4X,3HCM5,4X,3HCM),4X,  3HCM5 , 4X  ,  3H(  M  5  ) 

W ;  .  '  •  ,  7  5 

x'P  :  5,  Vjj  T  I  ME  (  1  5  ,  ANGPDC  1  5  ,  ANGRSDC 1  5  ,  ANGUEL  (15,  ACHC  1  5  ,  ACU(  1  5  , 


FORCE  Cl  )  ,  AChCGC  1 )  ,  ACUCGC  1  )  ,  ACXHUC  1  )  ,  ACYHUC  1 )  ,  XACRSC  1  )  , 

XACPOT  C 1 5 , XRPOT Cl), YACRSC 1 ) , YACPOT  Cl), YRPOT C 1 ) 

FORMAT  i FB . 5 , 5F8 . 3 , FB . 2 , 4FB . 3 , 6F7 . 3 ) 

WRITE  (6,15  1  TIME  C I  +  1 ) , ANGPDC I  +  1 ) , ANGRSDC2 ) , ANGUELC I  +  1 ) , ACHC I +1 ) , 
ACUC 1  +  1 ) , FORCE  C 1  +  1 ) , ACHCGC 2) , ACUCGC 2) , ACXHUC 2) , ACYHUC2) , 
XACRSC2)  ,  XACPOT  (8)  ,  XRPOT C  1  +1  )  ,  YACRSC 2 )  ,  YACPOTC 2 )  ,  YRPOTC  I  +  1  ) 

SET  COMPUTED  UALUES  AS  INITIAL  UALUES  FOR  THE  NEXT  ITERATION 

ANGRSnCl)  =  ANGRSDC2) 

ANGRS Cl'1  -  ANGRSC2 ) 

ACHCGC1)  -  ACHCGC2) 

ACUCGC 1)  -  ACUC3C2) 

ACXHUC1)  -  ACXHUC2) 

ACYHUC1)  -  ACYHUC2 ) 

UELACXC1)  -  UELACXC2) 

UELACYC1)  -  UELACYC2) 

XACRSC1)  -  XACRSC2) 

YACRSC 1)  -  YACRSC2) 

ANGPOTC1)  “  ANGP0TC2) 

ACHCGPC 1 )  -  ACHCGPC2) 

ACUCGPC 1 )  -  ACUCGPC2) 

ACXHUPC1)  -  ACXHUPC 2 ) 

ACYKUPC1)  -  ACYHUPC2) 

UELAXPC 1 )  -  UELAXPC2) 

UELAYPC1)  -  UELAYPC2) 

XACPOTC1)  -  XACP0TC2) 

YACPOTC 1)  -  YACP0TC2) 

IF  CI.EQ.NSTOP)  STOP 
I  -  I  +  l 


LIST  OF  SYMBOLS 


acceleration  (underlining  denotes  a  vector  qunntitv) 

acceleration  output  of  the  initially  horizontal  accelerometer 

acceleration  due  to  rotation  sensed  by  initial Iv  horizontal 
accelerometer 

starting  acceleration  for  the  horizontal  accelerometer 
aj| ]  corrected  for  gravity 

acceleration  at  time  t2  for  the  horizontal  accelerometer 
a^2  corrected  for  gravity 

acceleration  magnitude  due  to  rotation  at  radius  R2  that  is 

along  radius  R2  and  normal  to  the  path  of  movement 

acceleration  magnitude  due  to  rotation  at  radius  R2  that  is 

perpendicular  to  the  radius  and  tangent  to  the  path  of  movement 

acceleration  output  of  initially  vertical  accelerometer 

acceleration  due  to  rotation  sensed  by  initially  vertical 
accelerometer 

starting  acceleration  for  the  vertical  accelerometer 
ay  1  corrected  for  gravity 

acceleration  at  time  t2  for  the  vertical  accelerometer 
ay2  corrected  for  gravity 

starting  horizontal  (X)  component  of  acceleration  at 
accelerometers  in  fixed  coordinate  system 

horizontal  (X)  component  of  acceleration  at  accelerometer  at 
time  t;  in  fixed  coordinate  system 

starting  vertical  (Y)  component  of  acceleration  at 
accelerometers  in  fixed  coordinate  system 

vertical  (Y)  component  of  acceleration  at  accelerometer  at  time 
ta  in  fixed  coordinate  system 

acceleration  due  to  gravity 

constant  length  of  force  moment  arm 

moment  of  inertia  about  axis  of  quarter  wheel  and  sensor  . 


-’I  '  .  i: c  mibtdnt  l  or  high  force  spring 
r  i  c-'iisr  n’t  for  low  force  spring 
"  :  quarter  wheel  and  sensors 

’  c  re  vv 1  ut  ions  per  minute  of  spin  table 
i :en  or  point  P  as  measured  in  the  body-centered  system 
;  a:  t  ■  t  he  seismic  mass  in  inches  on  a  spin  table 

;es.;  L'  n  ot  point  P  in  fixed  coordinate  system 
initial  position  of  point  P  in  fixed  coordinate  system 
arci  'oration  of  point  P  in  fixed  coordinate  system 

r-.ci’is  from  axis  to  sensitive  masses  of  the  accelerometers  on 

t’ne  quarter  wheel 

position  in  fixed  coordinate  system  of  origin  of  body-centered 

cooi d inate  system 

-  radius  from  axis  to  center  of  gravity  of  quarter  wheel 
radius  from  axis  to  accelerometers 

.tree  1  er.it ion  of  the  origin  of  the  body-centered  system 

c  ;  to. 

-  vr  time 

c a  odd  rime  value 

■■  component  of  velocity  at  accelerometers  at  start  time 

X  component  of  velocity  at  accelerometers  at  time  t£ 

-  .  •  ■mponent  of  velocity  at  accelerometers  at  start  time 

■  Ticnt  of  velocity  at  accelerometers  at  time  t2 

i-  if.  of  point  P  in  a  rigid  body  relative  to  fixed  coordinate 

.  --r.f'n 

; •»  1  velocity  of  point  P  in  fixed  coordinate  system 
i  t.  ,  *ii  the  quarter  wheel  and  sensors 

■■'ing  horizontal  displacement  of  the  accelerometers  in  fixed 

.  i  •  1  i  n.i  t  e  s  v  s  t  em 
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X2  =  horizontal  displacement  of  the  accelerometers  at  time  t2  in 
fixed  coordinate  system 

Y i  =  starting  vertical  displacement  of  the  accelerometers  in  fixed 
coordinate  system 

'i 2  ~  vertical  displacement  of  the  accelerometers  at  time  t2  in  fixed 

coordinate  system 

Z|  =  initial  extension  in  length  beyond  its  rest  position  of  high 
force  spring  in  test  facility  model 

7-2  =  initial  extension  in  length  beyond  its  rest  position  of  low 

force  spring  in  test  facility  model 

0  =  angle  of  rotation  from  an  initial  position 

0 o  =  initial  angle  of  R  with  respect  to  the  horizontal 

0j  =  initial  value  of  rotation  angle  0  (0^  =0),  or  initial  angle  of 

Rj  with  respect  to  the  horizontal 

02  =  initial  angle  of  R2  with  respect  to  the  horizontal 

0rS 1  =  initial  angle  of  rotation  of  rate  sensor 

0RS2  =  integrated  value  for  angle  of  rotation  at  time  t2  using  angular 
velocity  from  rate  sensor 

0p2  =  rotation  angle  at  time  t2  derived  from  potentiometer  record 

0  =  angular  velocity 

®RS1  =  initial  angular  velocity  from  rate  sensor 

®RS2  ~  angular  velocity  at  time  t2  from  rate  sensor 

d^9 

9  =  — angular  acceleration 

dt" 

w  =  angular  velocity  of  the  body-centered  system  relative  to  the 

fixed  system 

uj  =  time  rate  of  change  of  angular  velocity 
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